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EXECUTIVE  SUMMARY 


This  report  summarizes  the  results  of  an  ongoing  SRI  International 
research  program  directed  toward  understanding  the  processes  that  control 
the  evolution  of  plasma  structure  in  both  the  natural  and  the  nuclear 
environments.  The  major  emphasis  during  the  past  year  has  been  on  high- 
latitude  processes,  because  the  high-latitude  region  is  poorly  understood 
compared  with  the  equatorial  ionosphere  where  the  conditions  and  occur- 
rance  of  structure  development  are  now  well  established.  Nevertheless, 
some  of  the  theoretical  developments  regarding  E-  and  F-region  coupling, 
which  were  initiated  because  of  the  importance  of  the  auroral  E  region, 
have  been  found  to  have  applications  to  equatorial  as  well  as  barium 
cloud  phenomena.  In  particular,  we  have  shown  theoretically  that  the 
formation  of  low-altitude  image  striations  is  highly  scale- size  selective. 
When  the  image  formation  process  is  driven  by  equatorial  spread-F,  the 
image  growth  rate  peaks  at  approximately  1  km.  We  have  shown  that  this 
is  in  good  agreement  with  rocket  observations  of  image  striations  at  the 
equator.  Further,  these  results  support  the  speculation  based  on  NRL 
simulations  that  E-region  image  striations  play  a  role  in  the  "freezing" 
of  barium  clouds.  Such  phenomena  are  important  for  determining  striation 
lifetime  and  size  distribution. 

Another  important  E-  to  F-region  coupling  effect  at  the  equator  is 
E-region  control  of  the  occurrence  of  equatorial  spread  F.  Following  a 
suggestion  of  Major  Leon  Wittwer  and  Dr.  Sidney  Ossakow,  we  have  shown 
that  the  seasonal  peaks  in  equatorial  spread-F  activity  occur  during 
months  when  sunset  is  nearly  simultaneous  at  the  conjugate  E  layers. 

During  these  periods  conductivity  gradients  are  changing  most  rapidly. 

Because  nonlinear  plasma  phenomena  couple  structure  at  different 
scale  sizes,  the  entire  spectrum  of  plasma  structure,  ranging  from  tens 
of  kilometers  to  centimeters,  must  be  studied  as  a  whole.  This  compli¬ 
cates  the  experimental  picture  because  one  must  combine  observations 
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from  separate  experimental  techniques,  each  of  which  can  survey  only  a 
limited  range  of  scale  sizes.  Thus,  at  large  scales  (wavelengths  >  10  km), 
we  have  shown  that  the  majority  of  F-layer  electron  density  enhancements 
("blobs”)  have  limited  both  zonal  and  meridional  extent.  The  structure 
observed  in  detail  for  one  such  blob  was  consistent  with  that  expected 
from  the  gradient  drift  instability.  In  the  intermediate-scale  regime 
(wavelength  <  1  km)  we  have  shown  that  the  irregularity  anisotropy,  both 

t  • 

perpendicular  and  parallel  to  the  magnetic  field,  has  a  strong  latitude 
and  local-time  variation.  The  anisotropy  appears  to  be  controlled  by  the 
magnetospheric  convection  pattern. 

An  ingiortant  experimental  breakthrough  made  during  the  past  year  is 
a  refinement  of  the  spaced- receiver  technique  to  determine  simultaneously 
irregularity  anisotropy  and  drift.  This  technique  provides  an  inexpensive 
means  of  continuously  monitoring  plasma  convection  at  high  latitudes. 
Moreover,  in  conjunction  with  HILAT,  it  provides  an  ideal  data  base  for 
refining  our  theoretical  models  of  high-latitude  irregularity  morphology 
developed  under  DNA  support. 
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I  INTRODUCTION 


The  motivation  for  studying  naturally  occurring  F-layer  irregular¬ 
ities  is  two-fold.  First  and  foremost,  these  irregularities  are  measur¬ 
able  and  from  their  measurements  we  can  extrapolate  to  the  more  severe 
propagation  disturbances  that  will  occur  in  a  nuclear  environment.  In 
this  regard,  naturally  occurring  F-layer  plasma  density  gradients  behave 
like  barium  clouds,  which  have  been  traditionally  used  to  simulate 
nuclear-produced  plasma  structure.  The  naturally  formed  plasma  density 
gradients  are  produced  at  high  latitudes  by  spatially  and  temporally 
structured  particle  precipitation  that  produces  local  plasma  density  en¬ 
hancements  or  "blobs"  in  the  F-region.  At  equatorial  latitudes,  gradients 
are  formed  on  the  edges  of  plasma  depletions  of  "bubbles"  formed  through 
the  Rayleigh-Taylor  process. 

In  the  case  of  the  blob,  or  bubble,  or  barium  cloud,  the  naturally 
occurring  driver  is  small  compared  with  the  driver  in  the  nuclear  case. 
However,  because  of  lower  plasma  temperatures  in  the  naturally  occurring 
driver,  the  diffusion  rate  is  also  much  smaller.  Therefore,  the  struc¬ 
ture  produced  in  barium  clouds  and  blobs  evolves  in  much  the  same  way  as 
the  nuclear-produced  structure.  Thus,  barium,  equatorial  bubbles,  and 
high-latitude  blobs  are  our  outdoor  laboratory  test  beds  for  models  of 
nuclear  plasma  structure  evolution. 

The  major  differences  among  these  models  lie  in  the  complexity  of 
the  plasma  systems.  The  barium  cloud  is  the  simplest  system  to  analyze. 
The  input  scale-size  is  known,  the  driving  forces  are  easily  determined, 
and  the  background  ionosphere  is  reasonably  well  understood.  The  equa¬ 
torial  ionosphere  represents  the  next  stage  of  conplexity.  The  system 
may  have  a  wide  range  of  scale  sizes  in  different  stages  of  evolution, 
as  determined  by  nature.  The  drivers  of  structure  in  the  F  layer  are 
neutral  atmospheric  winds,  gravity  waves,  and  the  dynamo  electric  field. 

At  high  latitudes  these  driv<"g  processes  are  accompanied  by  destabilizing 
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current  systems,  magnetospheric  convection  electric  fields,  and  dynamic 
background  ionospheric  E  and  F  layers.  If  a  nuclear  effects  code  were 
sufficiently  robust  to  reproduce  the  structure  observed  in  the  high- 
latitude  environment,  it  would  lend  confidence  to  its  ability  to  predict 
accurately  the  structure  resulting  from  a  nuclear  detonation. 

In  addition  to  code  verification,  the  second  motivation  for  studying 
naturally  occurring  plasma  density  irregularities  is  that  they  can  hinder 
DoD  communications  and  surveillance  systems.  Therefore,  understanding 
the  processes  that  control  the  evolution  and  ultimate  morphology  of 
naturally  structured  plasma  can  help  to  predict  system  performance  and 
to  develop  interference  mitigators. 

This  report  summarizes  the  results  of  an  ongoing  SRI  International 
research'  program  aimed  at  understanding  the  processes  that  control  the 
evolution  of  plasma  structure  in  the  natural  and  the  nuclear  environments. 
The  program  encompasses  both  theory  and  experiment  at  high  and  low 
latitudes.  Moreover,  because  a  vast  range  of  scale  sizes  (tens  of  kilo¬ 
meters  to  centimeters)  must  be  studied  simultaneously  to  a ssess  structure 
evolution  adequately,  a  rich  variety  of  complementary  experimental  tech¬ 
niques  must  be  employed  jointly.  Here  we  will  merely  outline  the  salient 
results;  the  detailed  experimental  procedures  and  theory  are  contained  in 
topical  reports  and  papers  to  which  we  shall  refer.  Section  II  describes 
irregularities  at  polar  latitudes,  and  Section  III  discusses  equatorial 
phenomena.  Section  IV  outlines  planned  future  work. 
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II  HIGH-LATITUDE  F-REGION  IRREGULARITIES 


An  important  finding  from  the  Wideband  satellite  experiment  was  the 
nighttime  scintillation  enhancement  that  occurs  with  high  probability 
whenever  the  propagation  path  lies  within  the  plane  of  the  local  L  shell 
[Fremouw  et  al..  1977;  Rino  et  al . .  1978}.  The  phenomenon  has  been 
attributed  to  a  geometrical  enhancement  from  sheetlike  irregularities. 
Evidence  is  accumulating,  however,  that  a  localized  source  region  also 
contributes  to  the  nighttime  scintillation  enhancement.  Through  combined 
incoherent-scatter  and  scintillation  observations,  the  F-region  source  of 
the  irregularities  and  the  dynamics  of  the  source  region  have  been  iden¬ 
tified  [Vickrey  et  al..  1980].  Enhanced  scintillation  appears  to  be 
associated  with  the  horizontal  plasma-density  gradients  ("walls")  of 
localized  plasma-density  enhancement  (blobs),  in  the  auroral  F  layer. 

The  experimental  and  theoretical  work  that  has  led  to  our  current 
understanding  of  this  high- latitude  phenomenon  is  reviewed  in  Rino  and 
Vickrev  [1982]  and  in  a  DNA  La  Jolla  meeting  report  to  be  published. 
Because  of  the  exceptionally  good  definition  of  the  source  region  and  its 
dynamics,  we  have  concentrated  on  this  particular  source  of  high-latitude 
irregularities.  The  results  are  grouped  into  two  scale  size  regimes 
according  to  the  following  criteria: 

•  The  structure  and  anisotropy  of  intermediate-scale 
(~  500  m  to  ~  10  km)  irregularities  that  are  the 
primary  contributors  to  the  scintillation. 

•  The  structure,  morphology,  and  likely  source  mechanisms 
of  large-scale  (>  10  km)  F-region  density  enhancements. 

The  physical  processes  that  couple  the  large-  and  small-scale  irre¬ 
gularities  are  discussed  in  detail  in  Kellev  et  al.  [1982]  and  Vickrev 
and  Kellev  [ 1982a],  and  have  been  reviewed  by  Vickrev  [1982]  and  Vickrev 
and  Kellev  [ 1982b].  These  authors  have  presented  a  framework  for  under¬ 
standing  high-latitude  irregularities  that  consists  of  (1)  a  simple  source 
function  of  irregularities,  (2)  a  model  of  the  scale-size  dependent  life- 
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time  of  irregularities  that  depends  on  E-region  conductivity,  and  (3)  a 
model  of  high- latitude  convection  that  redistributes  irregularities 
throughout  the  polar  regions.  Although  this  framework  is  straightforward, 
it  has  led  to  more  extensive  experiments  that  have  revealed  the  cooqplexity 
of  each  element.  In  the  following  two  sections,  we  highlight  the  perti¬ 
nent  results  of  some  of  those  experiments. 

* 

A.  Intermediate-Scale  Auroral  F-Regloo  Irregularities 

During  the  active  measurement  phase  of  the  Wideband  satellite  experi¬ 
ment  (May  1976  to  August  1979),  spaced-receiver  measurements  were  made  to 
resolve  unambiguously  the  anisotropy  and  drift  of  the  irregularities.  An 
elaborate  data  analysis  procedure  must  be  used  for  such  measurements 
[Rino  and  Livingston.  1982]  but  the  results  that  are  now  emerging  seem  to 
justify  the  effort.  However,  before  presenting  some  of  these  results,  we 
emphasize  that  the  correlation  method  is  most  sensitive  to  kilometer- scale 
structures.  Ultimately,  we  would  like  to  know  the  scale-size  or  spatial 
wavelength  dependence  of  the  anisotropy  and  drifts. 

Livingston  et  al.  [1982]  have  carefully  measured  the  variation  of 
the  axial  ratio  and  orientation  angle  of  the  diffraction  pattern  during  a 
large  number  of  satellite  passes.  They  were  able  to  identify  systematic 
transitions  among  three  generic  types  of  irregularity  structures:  rods, 
winds,  and  sheets.  Rods  are  isotropic  in  the  plane  transverse  to  the 
geomagnetic  field  and  can  be  characterized  by  an  axial  ratio  (aprallel  to 
transverse  dimensions)  of  a:l.  Wings  and  sheets  are  elongated  both  along 
the  geomagnetic  field  and  in  the  transverse  plane  along  the  direction  of 
the  local  L  shell  and  are  characterized  by  two  axial  ratio  parameters  as 
a:b:l.  For  sheets,  a«*b.  Wings  are  intermediate-scale  structures  with 
b  -  a/2. 

Figure  1  shows  the  general  pattern  established  by  Livingston  et  al. 
[1982].  The  contours  indicate  rods,  and  the  cross-hatched  regions  denote 
sheets  or  wings;  the  values  are  the  extension  ratio,  a.  The  axial  ratio 
of  the  rods  systematically  decreases  toward  the  poleward  boundary  of  the 
scintillation  zone,  giving  way  to  surprisingly  low  (nearly  isotropic) 
values.  Sheets  and  wings  are  confined  to  the  equatorward  region  of  the 
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FIGURE  1 


MORPHOLOGY  OF  IRREGULAR  ANISOTROPY  UNDER  CONDITIONS  OF 
MODERATE  ACTIVITY.  The  expected  large-scale  convection  pattern  is  superimposed. 


scintillation  zone.  There  is  a  systematic  reduction  in  anisotropy  (the 
rods  shorten  while  sheets  give  way  to  wings)  through  magnetic  midnight. 

Livingston  et  al.  [1982]  showed  that  this  pattern  is  closely  related 
to  the  expected  nighttime  convection  flow,  as  Figure  2  depicts.  The 
general  reduction  in  anisotropy  coincides  with  the  expected  duration  of 
the  Harang  discontinuity.  Sheets  occur  only  in  the  regions  of  strong 
electrojet  flows,  which  suggests  that  convection  and  velocity  shear  are 

r 

important  factors  in  controlling  the  irregularity  anisotropy  [Vickrey 
and  Kellev.  1982]. 

More  recent  measurements  using  polar  beacon  satellites  in  highly 
eccentric  orbits  provide  additional  evidence  for  the  relation  between 
convection  and  irregularity  anisotropy  (see  Section  Ill).  For  these 
satellites,  the  orbital  velocity  component  is  primarily  along  the  propa¬ 
gation  line  of  sight  and,  as  such,  does  not  strongly  affect  the  corre¬ 
lation  measurements. 

It  is  therefore  possible  to  establish  irregularity  drift,  as  well  as 
anisotropy  at  a  fixed  penetration  point.  An  example  of  anisotropy  varia¬ 
tion  for  a  few  hours  of  observation  is  shown  in  Figure  3  in  the  same  for¬ 
mat  as  Figure  1.  These  data  correspond  to  the  same  magnetic  conditions 
as  Figure  1,  and  the  local  time  variation  is  similar  to  that  in  the  long¬ 
term  morphological  pattern.  In  this  case,  however,  we  have  direct  con¬ 
firmation  from  the  velocity  data  in  Figure  4  that  the  anisotropy  minimum 
near  2300  corresponds  to  the  Harang  discontinuity. 

B.  Large-Scale  Auroral  F-Region  Irregularities 

To  complement  the  scintillation  observations  of  intermediate-scale 
auroral  structures  just  described,  the  Chatanika  incoherent-scatter  radar 
has  been  used  to  investigate  the  nature  of  the  large-scale  electron  den¬ 
sity  structure  found  in  the  auroral  F  layer.  The  radar  beam  was  scanned 
in  elevation,  usually  in  the  magnetic  meridian,  to  map  out  the  spatial 
characteristics  of  this  structure.  The  structure  has  been  found  to  take 
the  form  of  plasma-density  enhancements,  called  F-region  blobs.  We  have 
found  that  blobs  are  produced  by  two  possible  sources  of  soft-particle 
precipitation,  and  that  they  have  longitudinal  as  well  as  latitudinal 


structure. 
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Beside  specifying  the  polar  cusp  and  the  "soft  zone"  of  particle 
precipitation  [Eather,  1969]  as  source  regions  for  F-region  blobs  [e.g., 
Kellev  et  al. ,  1982] ,  we  have  found  that  at  least  two  specific  sources 
can  be  identified.  We  show  below  that  the  longest  (in  longitude)  blob 
is  located  near  the  equatorward  boundary  of  the  auroral  oval,  and  that 
the  shorter  blobs  are  located  poleward  of  this  boundary  blob.  Because 
of  its  characteristic  location  and  longitudinal  extent,  we  identify  the 
boundary  blob  as  a  unique' auroral  F-region  feature,  distinct  from  the 
other  ordinary  blobs.  We  suggest  that  the  source  of  the  boundary  blob 
is  latitudinally  localized  precipitation  of  low-energy  electrons  from 
the  inner  edge  of  the  central  plasma  sheet.  We  further  suggest  that  the 
source  of  the  other  blobs  is  the  inverted-V  event  [e.g.,  Frank  and 
Acker son.  1971]. 

We  present  in  Figure  5  examples  of  the  boundary  blob  in  a  set  of 
four  successive  scans  made  with  the  Chatanika  radar.  The  data  have  been 
plotted  as  isodensity  contours  in  altitude  and  ground  distance  coordinates. 
The  boundary  blob  in  Figure  5  is  the  plasma-density  (shaded)  region  that 
appears  to  move  equatorward  as  a  function  of  time.  Equatorward  of  the 
boundary  blob  is  a  sunlit  F  layer,  which  is  separated  from  the  boundary 
blob  by  a  narrow  latitudinal  region  of  low  plasma  density.  Poleward  of 
the  boundary  blob  are  the  ordinary  blobs.  Examples  of  isolated  blobs, 
like  those  in  the  third  panel  of  Figure  5,  suggest  that  the  blobs  are 
localized  enhancements  rather  than  depletions.  The  boundary  blob  is 
identified  by  its  proximity  to  the  equatorward  boundary  of  the  auroral 
oval,  defined  by  the  equatorward  edge  of  the  auroral  E  layer  that  is 
evident  in  all  four  of  the  panels  in  Figure  5. 

The  locations  of  the  boundary  blob  and  other  blobs  with  peak  plasma 
densities  greater  than  4  x  1GJ  el/cm^  have  been  plotted  as  a  function  of 
invariant  latitude  and  magnetic  local  time  (MLT)  in  Figure  6.  Solid 
circles  connected  by  line  segments  are  used  to  show  the  location  of  the 
boundary  blob.  The  longitudinal  extent  of  the  boundary  blob  is  evident 
from  its  persistence  as  a  function  of  MLT  and  from  comparison  with  the 
sporadic  appearance  and  disappearance  of  the  ordinary  blobs.  Electric- 
field  measurements  indicate  that  the  boundary  blob  is  longitudinally 
extended  and  is  not  an  isolated  feature  corotating  with  the  earth.  As  a 
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FIGURE  6  BOUNDARY  BLOB  LOCATION  AS  A  FUNCTION  OF  INVARIANT  LATITUDE 
AND  MAGNETIC  LOCAL  TIME 
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statistical  reference,  we  have  also  plotted  the  equatorward  boundary  of 
the  auroral  oval  (Q  ■  5),  as  determined  analytically  by  Starkov  [1969], 

We  see  from  Figure  6  that  the  invariant  latitude  of  the  boundary  varied 
as  a  function  of  MLT  in  a  manner  that  more  or  less  paralleled  that  of 
the  auroral  oval  boundary  in  the  afternoon  and  midnight  sectors.  The 
equatorward  departures  in  the  evening  and  postmidnight  sectors  appear  to 
be  associated  with  substQrm  effects. 

The  proposed  relationship  of  ordinary  blobs  to  inverted-V  events 
appears  reasonable  if  we  consider  the  nature  of  inverted-V  precipitation. 
Inverted-V  events  are  characterized  by  an  electron  energy  spectrum  whose 
peak  energy  varies  as  a  function  of  latitude  from  a  few  hundred  electron 
volts  to  a  few  kiloelectron  volts  before  decreasing  again  to  a  few  hundred 
electron  volts  le.g.,  Lin  and  Hoffman.  1979;  Hoffman  and  Lin.  1980]. 

Their  latitudinal  widths  vary  from  less  than  0.5  deg  to  no  more  than  2  deg, 
with  the  most  probable  width  occurring  at  the  smaller  values.  Inverted-V 
events  are  associated  with  intense  electron  fluxes  and  represent  the  most 
common  form  of  precipitation  in  the  auroral  zone.  For  example,  Hoffman 
and  Lin  [1981]  have  found  from  Atmospheric  Explorer-D  satellite  data  that 
inverted-V  events  occurred  in  94  percent  of  the  passes  analyzed,  with  72 
percent  of  those  events  occurring  near  the  poleward  edge  of  the  auroral 
oval.  The  spectral  characteristics  of  inverted-V  precipitation  whould 
produce  ionospheric  signatures  that  in  altitude-latitude  coordinates 
resemble  in  form  an  upright  V.  Because  of  more  rapid  loss  rates  at  lower 
altitudes,  the  upright-V  ionization  form  should  decay  into  a  pair  of  blobs. 
Of  course,  if  the  inverted-V  event  is  very  narrow  in  latitudinal  width, 
the  radar  signature  may  be  a  single  blob.  An  example  of  an  upright-V 
signature  is  presented  in  Figure  7.  The  isodensity  contours  indicate  the 
presence  of  an  auroral  arc  located  at  an  altitude  of  about  150  km.  If 
the  auroral  arc  were  produced  by  a  narrow  beam  of  keV  electrons,  we  would 
expect  the  plasma  density  to  be  latitudinally  confined  along  the  geomag¬ 
netic  field  line  passing  through  the  auroral  arc,  indicated  by  the  dashed 
line  in  Figure  7.  Instead,  we  see  a  plasma-density  contour  pattern  that 
spreads  outward  latitudinally  with  increasing  altitude,  as  is  expected 
from  an  inverted-V  event. 
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FIGURE  7  IONOSPHERIC  SIGNATURE  OF  AN  INVERTED-V  EVENT 


We  have  found  that  ordinary  blobs  can  have  longitudinal  structure. 

That  is,  blobs  are  not  only  limited  in  latitudinal  extent,  but  occasionally 
show  variations  in  latitudinal  cross  section  as  a  function  of  longitude. 
This  longitudinal  structure  is  consistent  with  plasma-density  structuring 
by  the  gradient-drift  instability  [e.g.,  Linson  and  Workman.  1970].  How¬ 
ever,  an  alternative  explanation  for  this  structure  is  a  varying  source 
of  soft-particle  precipitation.  This  second  source,  however,  appears  to 
require  very  large  particle  fluxes. 

A  Chatanika  data  set  used  to  reveal  the  longitudinal  structure 
associated  with  a  blob  is  presented  in  Figure  8.  The  data  set  consists 
of  six  elevation  scans  in  the  magnetic  meridian.  The  scan  sector  was 
centered  on  the  blob  and  limited  (compared  to  the  full  scans  shown  in 
Figure  5)  so  that  the  scans  were  closely  spaced  in  time.  The  isodensity 
contour  plots  in  Figure  8  differ  from  those  in  Figures  5  and  7  because 
the  data  were  reformatted  to  produce  a  plot  in  which  the  geomagnetic 
field  lines  are  vertically  oriented.  The  vertically  oriented  blobs  are 
therefore  field  aligned. 

Evidence  of  longitudinal  (or  east-west)  structure  is  contained  in 
the  variability  of  blob  characteristics  from  map  to  map.  From  Figure  8 
we  see  that  the  density  structure  found  in  successive  maps  does  not 
correlate  well.  A  blob  found  in  one  map  cannot  be  easily  identified  with 
a  blob  in  the  next  map  without  invoking  anomalously  large:  production  by 
particle  precipitation,  latitudinal  transport  by  east-west  electric  fields, 
or  loss  rate.  The  most  reasonable  interpretation  appears  to  be  in  terms 
of  east-west  structure.  (Details  of  this  data  set  are  described  in 
Tsunoda  and  Vickrey  [1982].) 

Interpretation  of  the  data  in  Figure  8  in  terms  of  longitudinal 
structure  is  best  illustrated  by  plan  views  of  the  isodensity  contours 
at  selected  altitudes.  To  make  these  plan  view’s,  we  have  taken  the 
plasma-density  variations  as  a  function  of  invariant  latitude,  at  selected 
altitudes  (horizontal  lines  in  the  maps  in  Figure  8),  and  fitted  isodensity 
contours  to  the  values  obtained  from  the  six  scans.  The  zonal  displace¬ 
ments  between  scans  were  computed  by  using  the  measured  southward  electric 
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field  of  22  mV/m.  No  attempt  was  made  to  compensate  for  any  latitudinal 
displacement  because  the  east-west  electric  field  component  was  small. 

The  resulting  isodensity  contour  maps  for  altitudes  of  360,  400,  and  440 
km  are  presented  in  Figure  9. 

The  apparently  uncorrelated  blob  features  seen  in  the  six  maps  in 
Figure  8  are  interpreted  in  Figure  9  as  longitudinal  structure  associated 
with  the  equatorward  wall  of  a  large-scale  blob.  Two  regions  of  enhanced 
plasma  density  have  "fingers"  extending  southward  and  are  separated  in 
east-west  direction  by  regions  of  depleted  plasma  density.  Similar 
patterns  of  east-west  structure  are  seen  in  the  three  contour  maps  taken 
at  different  altitudes.  The  east-west  separation  of  the  fingers  is  about 
150  km.  The  contour  maps  in  Figure  9  also  indicate  structuring  of  the 
west  wall  of  the  center  finger  in  Figure  9.  The  structuring  is  in  the 
form  of  an  eastward  penetration  of  a  low  plasma-density  contour  between 
two  regions  of  higher  plasma  density.  This  structure  is  most  evident  at 
the  400-km  altitude,  although  there  is  also  some  indication  of  the  east¬ 
ward  penetration  (of  a  lower  valued  contour)  at  the  360-km  altitude. 

The  blob  structure  described  above  is  consistent  with  that  produced 
by  the  gradient-drift  instability.  This  instability  has  a  positive 
growth  rate  whenever  the  plasma  drift  (relative  to  the  neutral  gas)  is 
parallel  to  the  plasma-density  gradient.  In  this  data  set,  the  electro¬ 
dynamic  conditions  included  a  southward  neutral  wind  of  370  m/s  and  a 
southward  electric  field  of  22  mV/m.  The  neutral  wind  acts  to  structure 
the  equatorward  wall  of  the  blob,  producing  the  longitudinal  structure 
just  described.  Similarly,  the  southward  electric  field  acts  to  struc¬ 
ture  the  west  wall  of  the  blob.  Both  forms  of  structure  were  observed 
as  shown  in  Figure  9. 

The  difficulty  with  this  interpretation  is  the  damping  effects  of  a 
conducting  E  layer.  If  the  E  layer  connected  to  the  blob  field  lines  was 
highly  conducting,  the  polarization  electric  fields  that  produce  the 
structuring  via  flux-tube  interchange  would  not  develop.  Unfortunately, 
the  E-layer  electron  density  was  not  measured  in  this  experiment  because 
the  blob  was  located  too  far  north  of  the  radar.  It  is  hoped  that  future 
experiments  will  resolve  this  ambiguity  in  interpretation. 
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The  magnetic  polar  cap  is  the  region  of  open  magnetic  field  lines 
that,  driven  by  the  solar  wind,  extends  into  the  tail  of  the  magneto* 
sphere.  It  is,  accordingly,  a  region  in  which  the  structure  and  dynamics 
are  closely  tied  to  the  state  of  the  interplanetary  magnetic  field  (IMF) 
and  thus  to  the  sun.  This  solar  control  is  conspicuous  if  we  consider 
the  scintillation  changes  observed  over  the  recent  sunspot  cycle.  A 
measurement  campaign  in  1975  judged  the  region  benign  in  terms  of  communi¬ 
cation  disruption,  but  in  1979  scintillation  levels  throughout  the  local 
winter  routinely  exceeded  those  simultaneously  observed  in  the  auroral 
zone  [Aarons  et  al..  1981]. 

The  polar  cap  region  is  interesting  not  only  because  of  the  severity 
of  scintillation,  but  it  is,  very  likely,  a  source  region  for  the  large- 
scale  electron  density  structures  observed  in  the  auroral  zone.  Further¬ 
more,  the  dynamics  of  the  region  provide  a  sensitive  indicator  of  the 
sense  and  magnitude  of  the  interplanetary  magnetic  field.  These  factors 
are  essential  to  an  understanding  of  the  irregularity  physics  and  to  the 
eventual  goal  of  global  irregularity  modeling. 

Until  the  observations  discussed  here  were  made,  there  were  virtually 
no  measurements  of  intermediate- scale  irregularity  anisotropy  in  the  polar 
cap.  The  occasional  exceptions  are  from  the  early  portions  of  Wideband 
satellite  passes  in  the  midnight  sector  and  under  disturbed  magnetic 
conditions;  hoover,  the  low»elevation-angle-observation  geometry  for 
these  Wideband  cases  is  far  from  ideal  for  the  detailed  diagnosis  of 
irregularity  shape.  In  contrast  with  this  paucity  of  anisotropic  data, 
there  has  been  extensive  study  of  plasma  drift  in  the  polar  cap  by  both 
modeling  and  observation.  Experimental  concentration  has  been  on  the 
large-scale  electric-field  patterns  that  are  mapped  down  from  the  magneto¬ 
sphere.  These  have  been  implied  from  multiple  magnetometer  data  [Friis- 
Christensen  et  al.,  1972]  or  directly  measured  by  balloon  flights  [Mozer 
and  Serlin,  1969]  and  satellites  [Heelis  et  al.,  1982], 


2 


The  polar  cap  plasma-flow  direction  is  dominated  by  the  northward 
(Bg)  component  of  the  interplanetary  field.  Under  normal  conditions, 
is  negative  and  small,  and  these  characteristics  cause  a  convection 
pattern  in  the  vicinity  of  the  magnetic  pole  that  is  two-celled,  with  flow 
from  noon  toward  midnight  and  sunward  return  through  the  dawn  and  dusk 
auroral  zones.  Generally,  the  two-flow  cells  are  asytanetric,  both  skewed 
from  the  noon-midnight  meridian  and  offset  from  the  magnetic  pole.  This 
lack  of  symmetry  is  generated  by  spatial  or  temporal  changes  in  the  east¬ 
ward  component  of  the  IMF,  By,  which  controls  the  dawn-dusk  flow  component 
to  first-order. 

When  Bz  turns  northward,  a  variety  of  complex  electric-field  distri¬ 
butions  can  occur.  Heppner  [1972]  has  categorized  these  in  terms  of  their 
dawn-dusk  signatures.  These  range  from  patterns  explainable  by  four- 
celled  polar-cap  convection  to  essentially  random  field  variations  through¬ 
out  the  entire  region.  Burke  et  al.  [1979]  have  studied  a  subset  of  Bz 
north  observations  and  have  postulated  a  magnetospheric  mechanism  for 
their  cause. 

Satellite  in-situ  measurements  of  polar-cap  electric  fields  like 
those  referenced  above  have  been  made  using  sun- synchronous  dawn-dusk 
orbits.  These  measurements  provide  one -dimensional  snapshots  of  the  most 
critical  features  of  the  convection  dynamics,  but  cannot  provide  simul¬ 
taneous  information  about  the  field  elsewhere  in  the  polar  cap.  The 
value  of  the  ground-based  drift  measurements  using  the  polar  beacons  is 
that  the  drift  pattern  can  be  continuously  monitored,  albeit  at  a  single 
point.  The  complementary  nature  of  the  two  methods  is  clear  and  should 
be  taken  full  advantage  of  during  the  HILAT  program. 

The  data  we  discuss  here  are  from  an  extensive  spaced -receiver  data 
set  collected  from  Thule,  Greenland.  Since  1979,  short  spans  of  multiple 
antenna  observations  have  been  routinely  made  daily  for  several  months  at 
a  time;  our  emphasis  here  is  on  data  from  several  days  of  continuous 
observations  made  during  January  1982.  The  diurnal  drift  patterns  derived 
from  these  data,  using  signals  from  the  UHF  polar  beacons,  are  currently 
being  consolidated.  Here  we  present  anisotropy  and  drift  data  for  one 
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consecutive  dawn,  noon,  and  dusk  sequence  during  which  the  magnetic-field 
conditions  were  nearly  constant. 

Dawn-sector  anisotropy  data  for  22  January  1962  are  shown  in  Figure 
10.  The  lower  scale  on  the  figure  is  universal  time;  the  upper  scale 
indicates  corrected  MLT.  Because  the  penetration  point  remains  nearly 
constant  during  the  observations,  the  two  time  scales  are  separated  by 
a  fixed  amount.  The  invariant  latitude  at  F-region  penetration  altitudes 
is  near  85  deg,  well  into  the  central  polar  cap. 

The  ground-diffraction-pattern  orientation  data  in  the  top  portion 
of  Figure  10  indicates  distinctly  rod-like  irregularities  as  the  cause  of 
the  scintillation.  In  fact,  the  model  predictions  cannot  be  distinguished 
from  the  measured  data;  those  for  sheet- like  irregularities  (not  shown) 
differ  by  nearly  30  deg  from  the  measured  orientation.  The  irregularity 
extension  for  which  two  model  predictions  are  shown  is  quite  variable, 
but  generally  is  greater  than  that  seen  in  the  auroral  zone  [Livingston 
et  al. .  1982],  reaching  20:1:1  beyond  0530  UT.  The  corresponding  dawn- 
sector  drifts  are  shown  in  Figure  11.  The  drift  magnitudes  are  high; 
they  exceed  900  m/s  slightly  past  0500  UT.  The  magnetic  north-south 
component  systematically  increases  from  200  m/ s  to  about  500  m/s  southward, 
while  the  more  rapidly  varying  east-west  component  averages  about  400  m/s 
in  an  eastward  direction.  For  this  invariant  latitude,  the  southeastward 
flow  is  consistent  with  an  overall  two-cell  flow  pattern  expected  for 
negative  conditions.  A  comparison  of  Figures  10  and  11  suggests  the  same 
convective  control  of  anisotropy  observed  in  the  auroral  zone  [Livingston 
et  al. ,  1982] .  Generally,  the  axial  ratio  follows  the  variations  in 
north-south  drift,  minimizing  when  the  drifts  are  high.  The  association 
between  anisotropy  and  the  east-west  velocity  component  is  not  clear. 

This  point  will  be  discussed  further  in  association  with  the  noon-  and 
dusk- sector  data. 

The  local-noon-sector  data  correspond  to  nearly  the  same  invariant 
latitude  of  85  deg  and  are  shown  in  Figures  12  and  13.  The  irregular¬ 
ities  are  still  rod-like,  but  are  extended  over  those  seen  in  the  dawn 
sector  and  exceed  50:1:1  at  times.  An  in-situ  extension  of  20:1:1  is  a 
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CORRECTED  MAGNETIC  LOCAL  TIME 
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FIGURE  10  LOCAL  DAWN-SECTOR  ANISOTROPY  DATA 
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CORRECTED  MAGNETIC  LOCAL  TIME 
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FIGURE  11  LOCAL  DAWN-SECTOR  DRIFTS 
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PATTERN 


FIGURE  12  LOCAL  NOON-SECTOR  DATA 
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typical  value.  The  axial  ratio  varies  with  a  60-  to  8°  ..an  period  up 
through  1300  UT,  perhaps  associated  with  the  passage  of  F-region  sun- 
aligned  arcs  [Weber  and  Buchau.  1981],  The  drifts  in  Figure  13  are  mostly 
antisunward,  as  expected  in  the  noon  sector,  but  also  include  a  distinct 
but  rapidly  varying  westward  component.  The  westward  component  presumably 
arises  from  an  asymmetry  in  the  two-celled  pattern  generated  by  a  nonzero 
component  of  By.  The  magnitude  of  the  drifts  occasionally  exceeds  800  m/s 
and  they  are  therefore  comparable  to  those  observed  in  the  dawn  sector. 
Unlike  the  dawn  sector,  however,  the  variations  in  axial  ratio  are  tied 
more  closely  to  changes  in  east-west  drift. 

Data  taken  near  local  dusk  on  the  same  day  are  shown  in  Figures  14 
and  15.  The  penetration  latitude  in  this  case  is  86  deg,  slightly  north¬ 
ward  of  the  other  time-sector  data.  The  irregularities  are  still  rod-like, 
systematically  increasing  from  about  20:1:1  elongation  in  midaftemoon  to 
about  40:1:1  by  2100  UT.  This  increase  is  associated  with  a  general 
decrease  in  drift  speed,  although  there  is  no  clear  association,  feature- 
for-feature,  with  the  northward  or  eastward  components.  The  northeast¬ 
ward  drift  in  this  case  also  fits  the  skewed  two-cell  flow  pattern  suggested 
by  the  dawn-aand  noon- sector  data. 

Detailed  comparisons  between  these  data  and  simultaneous  ionosonde 
and  optical  observations  are  under  way.  Associations  between  particular 
anisotropy  and  drift  features,  and  the  large-scale  polar-cap  structures 
can  then  be  established,  and  the  process  of  identifying  instability  pro¬ 
cesses  can  be  initiated.  From  these  measurements  alone,  however,  we  can 
make  some  preliminary  conclusions: 

(1)  The  irregularities  in  the  ~  600-m  spatial  wavelength 
regime  are  almost  exclusively  rod-like  and  range  in 
extension  from  about  5:1:1  to  50:1:1.  We  note  that 
on  the  day  following  the  examples  shown  here, 
deviations  from  the  rod- like  structures  do  occur 
that  correspond  to  unusual  drift  flow  patterns. 

For  either  irregularity  type,  scintillation  levels 
for  field-aligned  propagation  can  be  greatly  and 
consistently  enhanced.  For  the  polar  cap,  this 
will  be  at  a  near-vertical  geometry. 
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(2)  The  measured  drift  patterns  in  the  central  polar  cap 
fit  the  expected  two-celled  pattern,  and  can  exceed 
8u0  m/s  in  the  dawn,  noon,  and  dusk  sectors. 

(3)  There  is  a  consistent  association  between  drift  magni¬ 
tude  and  direction,  and  irregularity  anisotropy. 

Convection  dominance  of  anisotropy  has  been  clearly  demonstrated  in  the 
auroral  zone,  and  evidently  a  similar  process  is  ongoing  in  the  polar 
cap.  The  large-scale  dynamics  of  high-latitude  fields  are  relatively 
well  understood  and  are  effectively  modeled.  If  the  anisotropy/drift 
associations  can  be  defined,  their  definition  will  greatly  simplify 
global  prediction  of  irregularity  occurrence.  These  associations  will 
be  examined  directly  during  joint  experiments  using  the  Sondrestrom  rada 
HILAT,  and  continued  observations  of  these  polar  beacons. 

In  addition  to  these  relatively  detailed  observations  of  the  irregu 
larity  anisotropy  and  drift  and  of  their  interaction,  the  velocity  obser 
vatlons  can  be  combined  to  yield  samples  of  the  large-scale  convective 
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flow  patterns  in  the  central  polar  cap.  From  on?  day  to  the  next,  there 
are  significant  changes  in  the  observed  flow,  both  in  magnitude  and 
direction.  There  are  major  features  in  the  patterns  that  persist,  how¬ 
ever,  and  these  are  reflected  in  Figure  16,  in  which  six  days  of  contin¬ 
uous  drift  measurements  have  been  averaged. 

In  the  figure,  the  origin  of  each  vector  indicates  the  invariant 
latitude  and  corrected  MLT  of  the  F-region  ray  path  penetration.  As  can 
be  seen,  these  data  were  taken  well  into  the  central  polar  cap.  The 
satellite  orbits  change  only  slightly  from  day  to  day,  so  each  vector 
indicated  includes  a  minimum  of  ten,  and  usually  more  than  twenty,  indi¬ 
vidual  observations.  The  velocity  magnitudes  in  the  noon  sector  average 
500  to  600  m/s. 

There  are  several  features  to  note  in  the  pattern.  First,  the  flow 
is  generally  antisunward  in  the  noon  sector,  which  is  expected  for  the 
Bz  south  IMF  sense,  which  prevailed  during  these  six  days.  In  the  dawn 
and  predawn  sector,  the  flow  pattern  indicates  a  return  to  the  auroral 
zone,  even  at  this  high-  latitude.  This  behavior  suggests  that  the  polar 
cap  boundary  may  be  as  far  north  as  80  deg  invariant,  compared  with  the 
average  ~  75  deg  boundary  seen  in  the  summer  [Heppner.  1972] .  In  contrast 
to  the  flows  in  this  sector,  which  change  little  from  day  to  day,  the  dusk 
drifts  are  highly  variable.  Accordingly,  the  several-day  average  velocities 
are  reduced  in  magnitude.  The  most  interesting  feature  in  the  pattern, 
however,  is  in  the  premidnight  hours  near  the  magnetic  pole.  The  drift  in 
this  region  is  consistently  sunward  and  relatively  large. 

These  are  latitudes  and  times  that  are  not  probed  from  one  ground 
station  and  a  limited  number  of  satellites.  We  have  interpolated  the  flow 
into  these  regions  using  a  simple  flow  pattern,  as  indicated  by  the  arrows 
in  Figure  16.  The  overall  picture,  then,  is  one  with  antisunward  flow 
along  a  near  noon-to-midnight  meridian,  with  direct  return  to  the  dawn 
auroral  zone,  and  a  more  complex  return  to  the  evening  auroral  zone  in¬ 
volving  a  counterclockwise  whorl  near  the  magnetic  pole.  This  is  a  much 
more  complex  pattern  than  the  simple  two-celled  models  [Heelis  et  al . . 

1982]  of  convective  flow  would  suggest.  It  is,  however,  very  similar  to 


one  of  those  systematically  observed  by  the  OGO-6  satellite  during  the 
summer  season  [Heppner.  1972]. 

An  advantage  of  the  spaced- receiver  scintillation  measurements  of 
drift  is  that  the  irregularity  anisotropy  and  strength  are  simultaneously 
extracted.  We  are  currently  studying  the  morphology  of  these  two  para¬ 
meters  within  the  framework  provided  by  the  flow  pattern.  This  study 
will  provide  information 'crucial  to  the  definition  of  irregularity  source 
regions,  decay  rates,  and  structuring  mechanisms  in  the  polar  cap. 
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Ill  EQUATORIAL  IRREGULARITIES 


As  was  discussed  in  Section  II  B,  there  is  still  much  uncertainty 
about  the  role  o£  the  E-region  conductivity  in  determining  the  nature  of 
F-region  plasma  structure  in  the  auroral  zone.  On  the  one  hand,  a  highly 
conducting  E  layer  damps  interchange  instabilities  by  shorting  out  polar¬ 
ization  electric  fields;  on  the  other  hand,  the  presence  of  the  E  layer 
may  make  it  possible  to  tap  the  magnetospheric  free-energy  source  (that 
might  otherwise  be  unavailable)  via  field-aligned  currents.  Thus,  a 
structured,  conducting  E  layer  may  provide  an  additional  driver  for  F- 
region  instability  processes. 

Although  the  dynamics  of  the  low-latitude  E  region  are  significantly 
less  complicated  than  those  in  the  auroral  zone  and  are  better  understood, 
the  role  of  field-line  coupling  of  the  unstable  bottomside  equatorial  F 
region  to  the  conjugate  E  layer  below  is  also  poorly  understood.  In  the 
next  section  we  show  definitively  that  the  occurrence  of  E-layer  sunset 
controls  the  seasonal  and  longitudinal  behavior  of  spread-F  occurrence. 
Such  a  relationship  had  been  suggested  earlier  by  Major  Wittwer  [Private 
Communication,  1982],  As  in  the  auroral  case,  the  question  of  whether 
the  control  is  the  result  of  the  removal  of  a  strong  damping  process  or 
the  enhancement  of  a  driving  process  will  await  further  experiments. 

In  the  latter  part  of  this  section,  we  discuss  another  "nonlocal" 
phenomenon  in  which  the  E-region  plasma,  far  away  from  a  locally  unstable 
F-region  situation,  is  affected  by  the  F-region  instability.  We  present 
the  first  experimental  evidence  that  E-region  image  striations  form  in 
response  to  the  turbulent  electric  field  produced  in  the  unstable  F-layer 
above.  Our  interest  in  E-region  image  striations  was  prompted  by  the 
extensive  discussion  at  the  DNA  La  Jolla  conference  about  the  possible 
role  of  images  in  the  "freezing"  of  barium  clouds.  This  freezing  effect 
is  the  production  of  a  preferred  minimum  striation  scale  size.  Such 
freezing  might  be  expected  on  the  basis  of  numerical  simulations  [McDonald 
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et  al . .  1980].  However,  comparisons  of  barium  cloud  data  with  the  U- 
shaped  curve  [Ossakow  et  al . .  1981]  show  that  the  scale  size  predicted  by 
the  curve  (using  a  classic  ambipolar  diffusion  coefficient  and  the  cloud- 
to-background  Pedersen  conductivity  ratio,  m,  of  the  cloud)  is  too  small 
by  approximately  a  factor  of  10  to  30.  The  scale  size  predicted  by  the 
curve  is  given  by  L  =  Cm|(Dl/V^(m  +  l)V(m  -  1))  ],  where  DA  is  the  cross¬ 
field  diffusion  coefficient  and  V  is  the  slip  velocity  between  ions  and 
neutrals. 

To  make  the  curve  agree  with  barium  cloud  observations,  we  must  in¬ 
crease  the  value  of  m,  or  the  value  of  ,  or  both  values.  Ossakow  et  al. 
[1981]  have  suggested  that  m  is  increased  through  the  formation  of  E-region 
image  striations.  Alternatively,  some  form  of  anomalous  diffusion  may 
enhance  D^.  We  point  out  below  in  Section  III  B  that  the  image  formation 
process  itself  is  highly  scale-size  selective  and  produces  a  favored  scale 
size  of  ~  1  km;  i.e.,  very  close  to  the  observed  freezing  scale  of  barium 
clouds.  This  suggests  that  this  process  is  the  key  to  connecting  the 
freezing  scale  to  the  L-shaped  curve.  However,  comparison  of  more  detailed 
two-level  simulations  with  improved  observations  is  required  to  settle  the 
issue  of  whether  an  increase  in  m  alone  is  sufficient  to  explain  the 
freezing  effect  or  whether  anomalous  diffusion  is  also  important. 

A.  Control  of  Equatorial  Spread-F  by  E-Region  Conductivity 

The  occurrence  of  equatorial  scintillations  has  hitherto  revealed  a 
puzzling  pattern  that  appears  to  vary  with  both  season  and  longitude 
[Aarons.  1977;  Basu  and  Basu.  1981],  Scintillation  activity  has  been 
found  to  maximize  at  some  observing  stations  during  the  equinoctial  periods 
[e.g.,  Roster.  1972]  and  at  other  stations  during  the  local  summer  [e.g., 
Livingston.  1980] .  This  longitudinal  variability  in  seasonal  dependence  is 
not  strongly  controlled  by  solar  activity  [Chandra  and  Rastogi.  1972; 
Roster.  1972];  that  is,  the  year  of  observation  does  not  alter  the  basic 
seasonal  dependence  of  scintillation  activity  at  a  given  station  [e.g., 
Rastogi.  1980] .  Besides  differences  in  seasonal  behavior  among  observing 
stations,  satellite  and  multistation  observations  have  revealed  distinct 
longitudinal  variations  in  equatorial  irregularity  occurrence  during  a 
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given  season  [Sharma  and  Muldrew.  1975;  Basu  et  al. .  1976;  Aarons  et  al. . 
1980;  Muldrew.  1980]. 

Although  longitudinal  variation  in  scintillation  activity  may  be 
expected,  for  example,  because  of  geomagnetic  field  variations  with  longi¬ 
tude  [Burke  et  al..  1980;  Maruvama  and  Matuura.  1980;  Abdu  et  al. .  1981], 
a  local  seasonal  dependence  is  not.  Interpretation  in  terms  of  a  local 
seasonal  dependence  is  questionable  because  the  mechanisms  believed  res¬ 
ponsible  for  the  F-region  irregularities  that  produce  equatorial  scintil¬ 
lations  involve  entire  geomagnetic  flux  tubes  that  extend  into  both  con¬ 
jugate  hemispheres  [e.g.,  Tsunoda.  1980;  Weber  et  al. .  1983].  Aarons  et 
al.  [1980,  1981]  attempted  to  resolve  this  contradiction  by  hypothesizing 
that  asymmetric  hemispherical  variations  in  F-region  plasma  density 
(associated  with  the  equatorial  anomaly)  might  be  responsible  for  the 
apparent  local  seasonal  dependence.  That  is,  for  a  given  percentage 
fluctuation,  the  absolute  path-integrated  irregularity  strength  that  pro¬ 
duces  scintillations  would  be  stronger  for  a  higher  mean  plasma  density. 

Interpretation  of  the  seasonal  and  longitudinal  dependence  of  equa¬ 
torial  scintillations  should  be  possible  within  the  framework  of  existing 
theories  of  equatorial  irregularity  production.  The  mean  behavior  of 
scintillation  occurrence  (excluding  the  sporadic  day-to-day  occurrence 
pattern)  should  be  amenable  to  description  (and  prediction)  by  the  under¬ 
lying  physical  processes.  There  seems  little  doubt  that  the  occurrence 
of  equatorial  scintillations  is  controlled  by  the  same  large-scale  pro¬ 
cesses  that  control  the  occurrence  of  equatorial  spread-F  (ESF).  In 
other  words,  the  production  of  irregularities  responsible  for  both 
equatorial  scintillations  and  ESF  is  controlled  by  the  effectiveness  with 
which  polarization  electric  fields  can  develop  in  the  nighttime  equatorial 
F  layer.  The  irregularity  production  mechanism  is  believed  to  be  the 
generalized  gradient-drift  (or  Rayleigh-Taylor)  instability.  The  term 
"generalized"  is  used  to  incorporate,  within  a  single  instability  des¬ 
cription,  the  different  sources  of  free  energy  that  by  themselves  are 
often  referred  to  by  various  names.  Under  this  term,  we  include  the 
collisional  Rayleigh-Taylor  instability  driven  by  gravity  [e.g.,  Dungev. 
1956;  Haerendel.  197C,  nudson  and  Kennel.  1975;  Scannapieco  and  Ossakow. 
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1976],  the  current-convective  instability  driven  by  field-aligned  currents 
[e.g.,  Ossakow  and  Chaturvedi.  1979],  as  well  as  the  ordinary  gradient- 
drift  instability  driven  by  neutral  winds  and  electric  fields  [e.g., 

Linson  and  Workman.  1970;  Tsunoda.  1981;  Tsunoda  and  White.  1981;  Tsunoda 
et  al..  1982]. 

In  this  section  we  compare  the  occurrence  patterns  of  equatorial 
scintillations  reported  in  the  published  literature  with  the  times  of 
sunset  in  the  conjugate  E  layers,  and  then  show  that  both  the  seasonal 
and  longitudinal  dependences  of  scintillation  activity  are  controlled 
by  the  rate  of  change  in  the  integrated  E-region  Pedersen  conductivity. 
That  is,  the  occurrence  maxima  in  scintillation  activity  at  a  given 
longitude  are  shown  to  coincide  with  times  of  the  year  when  sunset  is 
simultaneous  at  the  conjugate  E  layers.  We  interpret  this  relationship 
as  support  for  the  concept  that  equatorial  irregularities  are  produced 
by  the  generalized  gradient-drift  instability,  driven  by  gravity  and  by 
field-aligned  currents  that  continue  to  flow  after  a  rapid  decay  of  the 
conjugate  E  layers,  as  well  as  by  an  eastward  neutral  wind. 


1.  Observations  of  Eouatorial  Scintillations 

Observations  of  equatorial  scintillations  reported  in  the  pub¬ 
lished  literature  are  summarized  in  Table  1.  All  scintillation  statistics 
except  those  presented  by  Taur  [1973],  have  been  reported  as  percentage 
occurrence.  Taur  [1973]  tabulated  the  average  monthly  number  of  minutes 
when  scintillations  exceeded  0.5  dB,  using  6-GHz  transmissions  from  vari¬ 
ous  INTELSAT  satellites.  Quantitative  comparisons  among  other  results 
(besides  those  of  Taur  [1973])  are  also  difficult  because  the  data  were 
obtained  at  different  frequencies  and  occurrence  statistics  were  computed 
using  different  references.  For  example,  Krishna  Moorthv  et  al.  [1978] 
used  scintillations  exceeding  0.25  dB  at  140  MHz,  Paulson  [1981]  used 
signal  fades  greater  than  6  dB  at  1.54  GHz,  and  Nichols  [1974]  used  signal 
fades  greater  than  5  dB  at  254  MHz.  The  Kwajalein  and  Ancon  results 
[Livingston.  1980]  are  based  on  the  scintillation  index  S4  exceeding  0.6 
at  137  MHz.  Other  factors  that  would  enter  into  a  quantitative  comparison 
are  the  elevation  angle  to  the  satellite  and  the  satellite  orbit.  Lower 
elevation  angles,  of  course,  would  produce  more  intense  scintillations. 
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OBSERVATIONS  OF  EQUATORIAL  SCINTILLATIONS 


Author( s) 

Place 

Frequency 

Observations 

Maxima 

Aarons  et  al. 

Accra,  Ghana 

257  MHz 

01/77  -  03/78 

Feb-Apr 

[1980] 

Huancayo,  Peru 

257  MHz 

01/77  -  03/78 

Feb-Apr 

Natal,  Brazil 

257  MHz 

01/77  -  03/78 

Nov-Jan 

Guam,'  U.S. 

257  MHz 

03/78  -  06/79 

Feb-Jul 

Bandvonadhvav  and 

Aarons  [1970] 

Huancayo,  Peru 

137  MHz 

07/67  -  02/69 

Sep -Nov 

Basu  et  al. 

[1980] 

Huancayo,  Peru 

1543  MHz 

04/76  -  10/77 

Mar, Oct 

Chatter iee  et  al. 

[ 1974] 

Huancayo,  Peru 

137  MHz 

05/67  -  05/68 

Oct-Nov, 

Mar 

Fane  and  Liu 

[1982] 

Hong  Kong,  U.K. 

4  GHz 

01/78  -  06/80 

Mar, Sep 

Koster  [1972] 

Legon,  Ghana 

137  MHz 

1967  -  1971 

Sep-Oct/ 
Mar- Apr 

Krishna  Moorthv 

et  al.  [1978] 

Trivandrum, 

India 

40,  140, 
360  MHz 

11/75  -  07/76 

Dec, Feb 
Jun 

Livineston  [1980] 

Ancon,  Peru 

VHF-UHF 

05/76  -  11/77 

Oct-Nov, 
Mar -Apr 

Kwajalein,  M.I. 

VHF-UHF 

10/76  -  10/77 

May/ 

Jul-Aug 

Nichols  [1974] 

Kwajalein,  M.I. 

254.1  MHz 

10/70  -  06/72 

May, Jul 

Guam,  U.S. 

258/1542 

MHz 

02/79  -  12/79 

Apr, Sep 

Taur  [1973] 

Hong  Kong,  U.K. 

4/6  GHz 

08/70  -  11/71 

Mar-Apr, 

Sep 

Guam,  U.S. 

4/6  GHz 

08/70  -  11/71 

Apr, Oct 

Tangua,  Brazil 

4/6  GHz 

08/70  -  11/71 

Feb-Mar 

Oct-Nov 

Walker  and  Chan 

[1970] 

Hong  Kong,  U.K. 

20,40  MHz 

10/64  -  07/68 

Summer 
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All  satellites  except  Wideband  [Livingston.  1980]  were  in  geosynchronous 
orbits.  Consequently,  scintillation  occurrence  maxima  could  be  determined 
from  continuous  measurements.  Wideband  was  in  a  (sun  synchronous)  noon- 
midnight  polar  orbit,  thus  limiting  scintillation  measurements  to  a  narrow 
local  time  sector.  Fortunately,  we  are  concerned  only  with  the  relative 
occurrence  frequency  of  scintillations,  on  a  monthly  basis  throughout  the 
year.  All  geosynchronous  satellite  results,  therefore,  can  be  compared 
directly;  those  from  the  Wideband  satellite  may  require  some  qualification. 

The  observations  can  be  categorized  into  two  longitudinal  regions 
the  American-African  sector  (which  includes  Huancayo  and  Ancon,  Accra  and 
Legon,  Tangua,  and  Natal)  and  the  Indian-Pacific  sector  (which  includes 
Trivandrum,  Hong  Kong,  Guam,  and  Kwajalein).  Scintillation  observations 
made  at  a  given  longitude  were  found  to  have  a  generally  repeatable 
seasonal  dependence;  e.g.,  observations  were  made  from  Huancayo  (including 
Ancon)  during  1967-69,  1976-77,  and  1977-78  showed  scintillation  activity 
maxima  from  February  to  April,  and  from  September  to  November.  The  repeat¬ 
ability  of  the  seasonal  dependence  in  scintillation  activity  indicates 
that  solar  activity  is  not  a  major  factor. 

In  the  following  section,  we  use  the  seasonal  dependences  found 
by  some  of  the  researchers  listed  above  to  establish  the  relationship 
between  scintillation  activity  and  the  sunset  times  at  the  conjugate  E 
layers.  A  difficulty  with  the  published  results  was  a  tendence  for 
researchers  to  reduce  scintillation  occurrence  data  into  three-month 
averages.  Averages  over  this  length  of  time  allow  interpretation  only  in 
terms  of  equinox  and  solstice,  and  tend  to  obscure  the  more  subtle  relation 
ship  between  the  solar  terminator  and  the  geomagnetic  flux  tubes  involved 
in  irregularity  generation.  Most  of  our  analysis  is  performed  using 
scintillation  activity  published  as  monthly  occurrence  frequency;  however, 
some  of  the  quarterly  averaged  results  are  used  in  cases  where  the  results 
from  that  station  were  considered  important  for  our  purposes. 

2.  Scintillation  Activity  and  E-Region  Sunset 

We  now  examine  the  seasonal  variation  in  sunset  times  at  the 
two  conjugate  E  layers  for  various  longitudes.  As  discussed  above,  we 


selected  four  locations  in  the  American-African  sector  (Huancayo,  Tangua, 
Nat»l,  Accra)  and  four  in  the  Indian-Pacific  sector  (Trivandrum,  Hong 
Kong,  Guam,  Kwajalein).  We  investigated  the  role  of  E-region  sunset  in 
the  seasonal  occurrence  of  equatorial  scintillations  be  determining  the 
local  solar  time  (LST)  at  the  dip  equator  when  each  of  the  conjugate  E 
layers  enters  darkness.  For  this  purpose  we  used  a  simplified  geometry 
in  which  the  geomagnetic  field  line  passing  over  a  given  observation 
station  is  characterized  by  the  geographic  latitude  at  which  it  crosses 
the  dip  equator  and  the  magnetic  declination  there.  The  conjugate  E 
layers  are  then  assumed  to  lie  at  fixed  distances  along  a  great  circle 
path  on  either  side  of  the  dip  equator.  This  geometry  is  presented  in 
Figure  17  in  earth-centered  spherical  coordinates.  The  locations  of  the 
north  geographic  pole  (N),  the  dip  equator  (Q),  the  subsolar  point  (S), 
the  conjugate  E  layers  (E,  E 7  ) ,  and  associated  angles  are  shown. 

We  are  interested  in  determining  the  local  solar  times  at  the 
dip  equator  when  the  solar  zenith  angle  at  each  conjugate  E  layer  becomes 
90  deg.  A  nominal  solar  zenith  angle  of  90  deg,  or  ground  sunset,  was 
selected  because  the  integrated  Pedersen  conductivity  in  the  E  layer 
approaches  the  nighttime  value  around  that  time  [Baron  et  al. .  1982; 
Senior.  1980;  Vickrey  et  al. .  1981;  Watt.  1973].  We  first  determine  the 
geographic  latitude  of  the  E  layer  by  using  the  law  of  cosines  on  triangle 
NEQ, 


0£  =  sin  ^(sin  9q  cos  l  +  cos  0^  sin  l  cos  D) 


(1) 


The  difference  in  hour  angle  (or  LST)  between  the  dip  equator 
and  the  E  layer  is  given  b  • 


6  55  sin”*(sin  l  sin  D/cos  0E) 


(2) 


Using  Eqs.  (1)  and  (2),  we  can  determine  the  solar  zenith  angle 
at  the  E  layer  from  the  following, 


X-.  *  cos  HsinQ  sin?,.  +  cos©  cos©_  cos(A+6)l 
t  st>s  t 


(  V 
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FIGURE  17 


EARTH-CENTERED  SPHERICAL  COORDINATES  RELATING  THE 
SOLAR  ZENITH  ANGLE  (x)  AT  THE  MAGNETIC  DIP  EQUATOR 
(Q)  AND  THE  CONJUGATE  E  LAYERS  (E,  E')  TO  THE  GEO¬ 
MAGNETIC  FIELD  GEOMETRY  (D,  8) 


When  the  solar  zenith  angle  at  the  E  layer  is  90  deg,  Eq.  (3) 


■.educes  to 


A  =  cos  1  (-tan8  tan9„)  -  6  (4) 

s  t. 

Equations  (1),  (2),  and  (4)  are  applied  to  both  E  layers  to 
determine  the  seasonal  dependence  of  LST  at  which  the  two  sunsets  occur. 

The  two  features  of  interest,  the  seasonal  variation  of  the 
times  when  the  conjugate  E  regions  enter  darkness  and  the  time  of  year 
when  both  E  regions  enter  darkness  together,  cannot  be  easily  inferred 
from  visual  inspection  of  the  above  equations.  We  therefore  applied  the 
above  equations  to  the  eight  selected  stations.  The  characteristics  of 
the  geomagnetic  field  lines  passing  over  the  eight  stations  are  listed 
in  Table  2. 


Table  2 

GEOGRAPHIC  AND  GEOMAGNETIC  PARAMETERS  OF  STATIONS 


Station 

Geog raphic 
( la t , long; 

Location 
in  deg) 

Dip  Equator 
( la t, long;  in  deg) 

Magnetic 

Declination 

(deg) 

Acc  ra ,  Ghana 

5.3, 

-1.1 

9.4, 

-0.7 

-7.8 

Guam,  U.S. 

12.4, 

147.0 

8.0, 

144.8 

3.2 

Hong  Kong,  U.K. 

22.2, 

114.2 

9.0, 

114.3 

0.5 

Huancayo,  Peru 

-11.0, 

-68.5 

-12.9, 

-75.1 

2.8 

Kwajalein,  M.I. 

9.4, 

167.5 

5.8, 

166.4 

8.8 

Natal,  Brazil 

-5.6, 

-33.7 

-2.8, 

-37.8 

-19.2 

Tangua,  Brazil 

-22.7, 

-42.8 

-8.0, 

-47.1 

-15.0 

Trivandrum,  India 

8.5, 

76.9 

9.5, 

76.7 

-2.9 

The  results  for  the  Indian-Pacific  sector  are  presented  in 
Figure  18,  The  family  of  oscillating  curves  in  each  panel  shows  the 
local  times  of  conjugate  E-region  sunsets  as  a  function  of  day  of  year. 
Each  curve  has  been  labeled  with  the  assumed  great-circle  ground-distance 
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FIGURE  18  COMPARISON  OF  THE  TIMES  OF  EQUATORIAL  SCINTILLATION  MAXIMUM 
TO  THOSE  OF  THE  SUNSET  NODES.  INDIAN-PACIFIC  SECTOR 
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separation  of  the  E  regions  from  the  magnetic  dip  equator.  Four  E-region 
spacings  (5,  10,  15,  20  deg)  were  selected  to  represent  expected  locations 
of  the  magnetic  flux  tubes  filled  with  scintillation-producing  irregular¬ 
ities.  All  the  curves  are  characterized  by  two  nodes  when  the  local  times 
of  E  region  sunset  are  identical  and  by  the  amplitude  of  their  local-time 
variation.  Of  particular  interest  are  the  days  of  the  year  when  the 
sunset  nodes  occur.  The  node  spacings  in  Figures  18(a)  through  18(d)  are 
seen  to  move  systematically  toward  the  June  solstice,  a  trend  produced  by 
the  increasing  (eastward)  magnetic  declination.  The  maximum  peak-to-peak 
amplitude  of  the  curves,  similar  for  all  four  stations,  is  primarily 
associated  with  the  geographic  latitude  of  the  conjugate  E  regions. 

The  occurrence  frequencies  of  equatorial  scintillations  observed 
at  each  of  the  four  stations  are  also  plotted  with  the  sunset  curves  in 
Figure  18.  Referring  first  to  Figure  18(b),  we  find  a  simple  relationship 
between  E-region  sunset  times  and  the  occurrence  maxima  for  scintillations 
i.e.,  scintillation  activity  appears  to  maximize  at  times  of  the  year  when 
the  sunset  nodes  occur.  Although  only  the  results  from  Taur  [1973]  are 
shown  in  Figure  18(b),  those  obtained  by  Fang  and  Liu  [1982]  are  nearly 
identical.  Similar  results  were  obtained  for  Guam,  as  shown  in  Figure 
18(c).  Although  the  6-GHz  results  [Taur.  19 73]  indicate  a  displacement 
between  a  scintillation  occurrence  maximum  and  the  fall  node,  the  1.5-GHz 
scintillation  measurements  obtained  by  Paulson  [1981]  indicate  coincidence 

The  Trivandrum  results  [Krishna  Moorthv  et  al..  1978]  shown  in 
Figure  18(a)  also  indicate  a  similar  relationship  between  the  maxima  in 
scintillation  activity  and  the  sunset  nodes.  The  scintillation  occurrence 
is  seen  to  peak  around  the  spring  node.  No  data  were  available,  however, 
to  indicate  whether  a  similar  maximum  is  associated  with  the  fall  node. 

To  supplement  the  scintillation  results,  we  have  plotted  curves  of  ESF 
occurrence  at  2100  LST  found  in  ionograms  taken  at  Kodaikanal  [Sastri  and 
Murthv.  1975] .  During  early  evening,  range-type  ESF  is  known  to  predom¬ 
inate  and  to  correlate  well  with  the  occurrence  of  scintillations  [e.g., 
Huang.  1970] .  The  correspondence  between  ESF  occurrence  maxima  with  the 
two  sunsc*  nodes  is  clear.  The  only  discrepancies  appear  as  a  December 
maximum  in  scintillation  activity  and  as  a  January  1968  maximum  in  ESF 
occurrence. 
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Only  the  Kwajalein  results  obtained  by  Livingston  [1980]  appear 
to  suggest  a  more  complex  relationship  between  s”.r.aet  nodes  and  maxima 
in  scintillation  occurrence.  The  results  from  Nichols  [1974]  indicate  a 
a  broad  maximum  centered  more  or  less  on  the  June  solstice.  The  broad 
maximum  suggests  that  conjugate  E-region  sunsets  do  not  have  to  occur 
simultaneously,  but  that  a  small  but  finite  time  difference  may  be  allowed 
The  structure  in  the  results  obtained  by  Livingston  [1980],  however,  seems 
to  suggest  a  close  but  not  coincident  relationship  between  scintillation 
occurrence  maxima  and  the  sunset  nodes.  The  primary  occurrence  maxima 
appear  to  be  slightly  displaced  toward  the  June  solstice;  there  is  even 
a  hint  of  local  minima  at  the  times  of  the  sunset  nodes.  Recall,  however, 
that  these  scintillation  measurements  were  made  in  the  midnight  sector  by 
a  polar-orbiting  satellite. 

The  relationship  between  seasonal  scintillation  activity  and 
conjugate  E-region  sunsets  for  the  four  stations  in  the  American-Af rican 
sector  are  presented  in  Figure  19.  The  Ancon  scintillation  results 
[Livingston.  1980],  shown  as  a  histogram  in  Figure  19(a),  reveal  a  coinci¬ 
dence  between  a  maximum  in  scintillation  activity  and  the  spring  node. 

The  other  scintillation  maximum,  however,  is  seen  to  be  shifted  away  from 
the  fall  node  toward  the  December  solstice.  Scintillation  results  from 
Huancayo  by  Chaturiee  et  al.  [1974]  are  seen  to  be  in  general  agreement 
with  the  results  from  the  polar-orbiting  Wideband  satellite.  The  third 
scintillation  results  in  Figure  19(a)  [Basu  et  al. .  1980]  also  support 
the  general  pattern  of  the  other  two  curves.  The  peaks  in  the  third 
curve,  however,  indicate  a  distinct  displacement  of  those  peaks  away  from 
both  nodes  toward  the  December  solstice.  There  appears  to  be  little  doubt 
therefore,  that  the  displacements  are  real  and  repeatable. 

The  other  results  presented  in  Figure  19  appear  to  be  in  accord 
with  our  finding  that  peaks  in  scintillation  activity  occur  during  the 
times  of  year  when  sunset  is  simultaneous  at  the  conjugate  E  layers.  The 
time  coincidence  is  striking  for  the  Tangua  results  [Figure  19(b)], 
although  there  is  a  slight  displacement  of  a  scintillation  maximum  away 
from  the  spring  node  towards  the  June  solstice.  The  comparisons  for  Natal 
and  Accra  also  reveal  this  relationship,  despite  the  fact  that  three-month 
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LOCAL  TIME  LOCAL  TIME 


(a)  HUANCAYO  (ANCON)  (b)  TANGUA 


DAY  OF  YEAR  DAY  OF  YEAR 

FIGURE  19  COMPARISON  OF  THE  TIMES  OF  EQUATORIAL  SCINTILLATION  MAXIMUM 
TO  THOSE  OF  THE  SUNSET  NODES,  AMERICAN-AFRICAN  SECTOR 
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averaged  scintillation  results  [Aarons  et  al..  1980]  were  used.  A  second 
scintillation  data  set,  obtained  by  Roster  L 1968] ,  also  reveals  the 
occurrence  of  scintillation  maxima  around  the  time  of  the  sunset  nodes. 

A  summary  of  the  relationship  between  scintillation  maxima  and 
sunset  nodes  is  presented  in  Figure  20.  We  have  plotted  the  months  of 
scintillation  maxima  and  the  sunset  nodes  for  each  of  the  eight  stations 
as  a  function  of  longitude:  Huancayo  (H),  Tangua  (T),  Natal  (N),  and 
Accra  (A),  in  the  American-African  sector;  and  Trivandrum  (Tr),  Hong  Kong 
(HK),  Guam  (G),  and  Kwajalein  (K),  in  the  Indian -Pacific  sector.  In 
addition,  we  have  supplemented  them  with  equatorial  scintillation  results 
[Taur.  1973]  from  Paumalu,  U.S.  ( P ) ;  Lurin,  Peru  (L);  Longovilo,  Chile 
(Lo);  Ascension  Island,  U.K.  (As);  Ras  Abu  Jarjur,  Bahrain  (R);  Djatiluhur, 
Indonesia  (D);  and  Taipei,  Republic  of  China  (Ta).(  At  the  Indian  longitude 
we  have  plotted  the  ESF  occurrence  results  from  Kodaikanal  [Sastri  and 
Murthv.  1975]  instead  of  the  scintillation  results  from  Trivandrum.  The 
longitudinal  variations  in  the  months  of  scintillation  maxima,  as  shown 
in  Figures  18  and  19,  are  seen  to  be  consistent  with  those  for  the  sunset 
nodes. 

The  monthly  trends  of  the  scintillation  maxima  and  sunset  nodes 
may  be  compared  with  the  longitudinal  variations  of  the  magnetic  declination 
and  geographic  latitude  of  the  dip  equator.  The  gentle  convergence  of  the 
scintillation  maxima  and  nodes  towards  the  June  solstice  in  the  Indian- 
Pacific  sector  is  seen  to  be  associated  with  a  corresponding  slow  increase 
in  positive  (eastward)  declination.  A  similar  variation  is  seen  in  the 
American-African  sector  where  the  scintillation  maxima  (and  nodes)  move 
toward  the  December  solstice  with  increasing  negative  (westward)  declina¬ 
tion.  The  effect  of  the  geographic  latitude  of  the  dip  equator  is  best 
seen  in  Figures  19(a)  and  19(c).  The  maximum  peak-to-peak  variation  in 
the  time  of  E-region  sunset  is  seen  to  be  smaller  for  stations  with  smaller 
separations  between  the  dip  and  geographic  equator  (Natal)  than  for 
stations  with  larger  separations  (Huancayo).  We  see  from  Figure  19(c) 
that  the  scintillation  maxima  merge  into  a  single  broader  peak  when  the 
nodes  are  closer  together  and  when  the  dip  equator  is  close  to  the  geo¬ 
graphic  equator. 
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FIGURE  20  COMPARISON  OF  THE  LONGITUDINAL  VARIATION  IN  THE  TIMES  OF  SCINTILLATION  ACTIVITY  MAXIMA  AND  SUNSET 
NODES  WITH  THE  MAGNETIC  DECLINATION  AND  GEOGRAPHIC  LATITUDE  OF  THE  MAGNETIC  DIP  EQUATOR 


Two  features  in  Figures  18  and  19  suggest  other  influencing 
factors  besides  simple  coincidence  in  sunset  times  of  the  conjugate  £ 
layers.  As  mentioned  earlier,  the  primary  scintillation  maxima  were 
occasionally  found  to  be  slightly  displaced  from  one  (or  both)  of  the 
sunset  nodes.  The  second  feature  is  an  asymmetry  in  the  occurrence  fre¬ 
quency  of  scintillations.  We  have  found  a  larger  peak  associated  with 
the  fall  node  than  with  the  spring  node  at  four  stations,  Hong  Kong,  Guam, 
Tangua,  and  Huancayo.  Paulson  [1981]  noted  the  asymmetry  in  his  Guam  data 
and  suggested  that  it  might  be  produced  by  varying  solar  activity.  We 
note,  however,  that  Taur  [1973]  also  obtained  a  similar  asymmetry  using 
data  taken  during  different  years.  A  third  feature,  the  isolated  occur¬ 
rence  of  a  scintillation  (and  ESF)  maximum  near  the  December  solstice  at 
Ancon  (and  Kodaikanal),  does  not  appear  to  be  real.  The  occurrence  fre¬ 
quency  at  Ancon  for  that  time  of  year  was  determined  from  only  two  satel¬ 
lite  passes  (the  average  number  of  passes  is  10)  [see  Livingston.  1980]. 

We  therefore  dismiss  that  estimate  as  spurious.  Not  enough  information 
is  available  on  the  Kodaikanal  results  to  assess  its  validity. 

Perhaps  the  most  revealing  piece  of  additional  information  is 
the  complete  absence  of  scintillation  activity  at  Huancayo  during  the 
June  solstice  when  darkness  at  both  E  regions  occurs  at  a  relatively 
early  local  time.  In  contrast,  there  is  significant  scintillation  acti¬ 
vity  at  other  times  of  the  year  when  the  last  sunset  at  a  conjugate  E 
layer  is  much  later  than  during  the  June  solstice.  This  finding  suggests 
that  simple  early  unloading  of  the  F-region  dynamo  [Rishbeth.  1971]  is 
not  sufficient  reason  for  enhanced  irregularity  generation  by  the  general¬ 
ized  gradient-drift  instability.  As  we  discuss  later  in  this  section, 
additional  (or  enhanced)  drivers  appear  to  be  present  during  the  sunset 
node  period  to  produce  a  significantly  increased  growth  rate  for  the 
generalized  gradient-drift  instability. 

3.  Sunset  Nodes  and  the  Longitudinal  Conductivity  Gradient 

Having  argued  that  simple  early  unloading  of  the  F-region  dynamo 
is  not  the  reason  for  enhanced  irregularity  generation,  we  now  suggest 
that  a  likely  source  of  free  energy  is  the  longitudinal  gradient  in  the 
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integrated  E-region  Pedersen  conductivity.  Several  researchers,  among 
them  Heelis  et  al.  [1974],  showed  that  this  conductivity  gradient  is 
capable  of  enhancing  the  eastward  electric  field  and  field-aligned  currents, 
both  of  which  are  drivers  for  the  generalized  gradient-drift  instability. 

We  show  that  longitudinal  gradients  are  steepest  around  the  times  of  maxima 
in  scintillation  activity. 

We  begin  by  deriving  formulas  for  the  longitudinal  gradient  and 
gradient  scale  length  of  the  integrated  E-region  Pedersen  conductivity. 
Senior  [1980]  showed  that  the  integrated  E-region  Pedersen  conductivity 
(produced  by  solar  radiation)  has  the  form 


Ep  ■  A  cosx  +  B 


(5) 


where  x  “  solar  zenith  angle,  A  =  12,  and  B  =  1  mhos  for  equinoctial 
periods.  We  assume  that  the  total  integrated  Pedersen  conductivity  is 
equal  to  the  sum  of  the  contributions  from  the  conjugate  E  layers.  Using 
Eqs.  (3)  and  (5),  the  total  integrated  E-region  Pedersen  conductivity  is 
given  by 

T 

Ep  *  2B  +  A  |  sinQ  s  (sin0£  +  sin0£/  ) 


+  cos0  [cos0„  cos(A+6)  +  cos0_  /  cos(A+6)]|  . 

St  t,  ' 

Differentiating  with  respect  to  longitude  yields 


-A  cos0  [cos0_  sin(A  +  6)  +  cos0_/  sin(A  +  6  /  )  ] 
st  t 


(6) 


(7) 


From  our  assumption  that  the  integrated  conductivities  are 
additive,  Eq.  (7)  shows  that  the  conductivity  gradients  are  also  addi¬ 
tive.  As  shown  below,  sunset  at  either  one  of  the  E  layers  therefore 
significantly  reduces  the  longitudinal  gradient  in  integrated  Pedersen 
conductivity. 
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length, 


Using  Eqs.  (6)  and  (7),  we  can  also  define  a  gradient  scale 


L 


(8) 


Although  a  gradient  scale  length  does  not  uniquely  characterize 
the  spatial  variation  in  conductivity,  we  find  it  useful  for  qualitative 
description  of  the  conditions. 

To  illustrate  the  behavior  of  the  longitudinal  conductivity 
gradient  [given  by  Eq.  (7)]  and  gradient  scale  length  as  a  function  of 
time  of  year,  we  computed  these  quantities  for  A  ■  90  deg,  i.e.,  sunset 
at  the  dip  equator.  We  further  assumed  that  the  integrated  Pedersen 
conductivity  decreased  according  to  Eq.  (5)  until  x  *  94.3  deg  when  Ep  = 
0.1  mho,  a  typical  nighttime  value  [e.g.,  Harper  and  Walker.  1977].  From 
this  time  on,  the  darkened  E  layer  is  assumed  to  contribute  no  longer  to 
the  longitudinal  gradient  or  gradient  scale  length. 

The  integrated  E-region  Pedersen  conductivity  gradient  and  the 
gradient  scale  length  for  each  of  the  eight  selected  stations  (Figures  19 
and  20)  are  plotted  in  Figures  21  and  22.  For  each  station  we  have  dis¬ 
played  the  gradient  scale  length  in  the  upper  panel  and  the  conductivity 
gradient  in  the  lower  panel.  As  before,  the  four  curves  correspond  to 
specified  spacings  of  the  conjugate  E  layers;  the  times  of  the  year  when 
the  sunset  nodes  occurred  are  indicated  by  vertical  lines.  The  node 
times  can  be  used  to  compare  these  curves  to  those  in  Figures  19  and  20. 

The  gradient  scale  length  curve  can  generally  be  described  by 
a  sinusoidal  variation  with  a  period  of  one  year,  interrupted  by  discon¬ 
tinuous  deviations  that  are  dependent  on  the  spacing  of  the  conjugate  E 
layers.  The  amplitude  of  the  one-year  cycle  is  proportional  to  the  geo¬ 
graphic  latitude  of  the  dip  equator.  For  example,  of  the  eight  stations, 
the  Huancayo  curve  has  the  largest  amplitude  and  the  Natal  curve  has  the 
smallest  amplitude.  The  phase  of  the  one-year  oscillation  is  dependent 
on  the  hemisphere  in  which  the  dip  equator  is  located:  those  in  the 
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INDIAN-PACIFIC  SECTOR 


(a)  TRIVANDRUM 


(b)  HONG  KONG 
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FIGURE  21  THE  VARIATION  WITH  TIME  OF  YEAR  IN  THE  LONGITUDINAL  GRADIENT 
AND  THE  GRADIENT  SCALE  LENGTH  ASSOCIATED  WITH  THE  INTEGRATED 
E-REGION  PEDERSEN  CONDUCTIVITY,  INDIAN-PACIFIC  SECTOR 
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Southern  Hemisphere  (Huancayo,  Natal,  Tangua)  have  minimum  gradient- 
scale  lengths  during  the  June  solstice;  those  in  the  Northern  Hemisphere 
(Accra,  Guam,  Hong  Kong,  Kwajalein,  Trivandrum)  have  minimum  during  the 
December  solstice. 

The  discontinuous  deviations  are  produced  by  the  abrupt  dis¬ 
appearance  of  contribution  to  the  gradient  scale  length  by  the  darkened 
E  layer.  When  this  happens,  the  gradient  scale  length  increases.  The 
increase  is  proportional  to  the  spacing  of  the  conjugate  E  layers  and 
inversely  related  to  the  geographic  latitude  of  the  dip  equator.  For 
example,  the  increase  in  gradient  scale  length  is  smallest  at  Huancayo 
and  largest  at  Natal.  Finally,  the  times  of  the  year  when  the  discontin¬ 
uities  occur  are  related  to  the  hemisphere  in  which  the  dip  equator  is 
located.  The  times  of  the  discontinuities  are  seen  to  move  toward  the 
December  solstice  for  stations  with  dip  equators  in  the  Southern  Hemi¬ 
sphere  (Huancayo,  Natal,  Tangua),  and  the  times  of  the  discontinuities 
move  toward  the  June  solstice  for  stations  with  dip  equators  in  the 
Northern  Hemisphere. 

The  curves  for  the  conductivity  gradient  are  simpler  in  charac¬ 
ter  than  those  for  the  gradient  scale  length.  They  are  characterized  by 
a  semiannual  oscillation  with  amplitude  peaks  during  equinoctial  periods 
and  discontinuous  decreases  produced  by  sunset  at  one  of  the  E  layers. 

The  semiannual  variation  in  conductivity  gradient  is  small.  The  seasonal 
variation  in  integrated  conductivity  gradient  therefore  appears  to  be 
controlled  by  the  occurrence  of  sunset  at  one  of  the  E  layers. 

The  peaks  in  scintillation  activity,  presented  in  Figures  19 
and  20,  are  found  to  occur  at  times  of  the  year  when  the  gradient  (gra¬ 
dient  scale  length)  in  integrated  E-region  Pedersen  conductivity  is 
steepest  (shortest).  This  conclusion  is  most  convincing  if  both  the 
gradient  and  the  gradient  scale  length  curves  are  used.  The  reader  is 
reminded  that  the  curves  in  Figures  21  and  22  refer  to  an  arbitrarily 
selected  time  (sunset  at  the  dip  equator)  and  are  intended  to  illustrate 
the  qualitative  behavior  of  the  longitudinal  gradient. 


From  comparisons  of  scintillation  activity  observed  at  the 
selected  eight  stations  with  sunset  times  at  the  conjugate  E  layers,  we 
have  reached  the  following  conclusions: 

(1)  Peaks  in  scintillation  (and  range-type  ESF)  activity 
occur  during  months  when  sunset  is  nearly  simultaneous 
at  the  conjugate  E  layers. 

(2)  Although  near- simultaneous  sunset  at  conjugate  E  layers 
might  also  be  interpreted  as  the  occurrence  of  earliest 
total  darkness,  scintillation  activity  is  not  related 
simply  to  months  of  early  darkness  at  both  conjugate 

E  layers.  The  Huancayo  results  provide  the  key  evidence 
for  this  subtle  differentiation. 

(3)  The  relationship  between  scintillation  occurrence  maxima 
and  sunset  nodes  is  not  exact.  Slight  displacements  (up 
to  about  a  month)  in  the  maxima  from  the  nodes  have  been 
found;  the  displacements,  amounting  to  less  than  10  per¬ 
cent,  are  not  necessarily  similar  for  the  spring  and 
f&xl  nodes. 

(A)  The  peak  value  of  scintillation  activity  associated  with 
the  fall  node  has  often  been  found  to  be  larger  than 
that  for  the  spring  node.  This  asymmetry  does  not  appear 
to  be  a  function  of  solar  activity,  as  suggested  by 
Paulson  [1981] . 

(5)  Peaks  in  scintillation  activity  coincide  with  the  times 
when  the  longitudinal  gradient  (gradient  scale  length)  in 
integrated  E-region  Pedersen  conductivity  is  steep  (short); 
conversely,  there  is  an  absence  of  scintillation  activity 
when  the  longitudinal  gradient  (gradient  scale  length)  is 
gradual  (long).  The  primary  variation  in  the  gradient 
value  is  produced  by  sunset  at  one  of  the  E  layers. 


B.  Low-Altitude  Image  Striations  Associated  with  Equatorial  Spread-F 

Ionospheric  plasma  instabilities  are  usually  discussed  in  terms  of 
local  parameters.  However,  as  was  illustrated  in  the  previous  section, 
because  electric  fields  of  scale  size  \  ^  1  km  map  along  magnetic  field 
lines,  plasma  populations  far  away  from  a  locally  unstable  region  can 
affect,  and  may  be  affected  by,  the  instability  process.  Another  example 
is  the  formation  of  "image"  striations  in  the  E  and  F^  regions  in  res¬ 
ponse  to  the  electric  field  produced  in  the  unstable  F  layer  above.  The 
image  results  principally  from  a  converging  velocity  field  operating  on 
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a  compressible  ion  gas.  In  this  section  the  theory  of  image  formation, 
including  * t:  scale-size  dependence,  is  discussed  in  detail.  Rocket 
measurements  of  electron  density  irregularities  are  also  presented.  We 
show  that  the  spectrum  of  density  structure  observed  off  the  magnetic 
equator  is  consistent  with  that  expected  from  the  image  formation  process 
(which  is  scale  size  and  density  dependent)  driven  by  the  gravitational 
Rayleigh-Taylor  instability  operating  on  the  bottomside  F  layer  at  the 
equator  itself. 

As  is  discussed  in  more  detail  below,  an  image  forms  at  altitudes 
where  the  ion  gas  is  compressible,  namely,  the  E  and  lower  F^  regions. 

To  date  such  features  have  not  been  verified  experimentally,  although 
computer  simulations  have  predicted  their  generation  [Goldman  et  al. . 

1976;  Lloyd  and  Haerendel.  1973]. 

The  data  discussed  here  come  from  fixed  bias  Langmuir  probes  flown 
on  two  rockets  launched  just  off  the  magnetic  equator  during  equatorial 
spread-F  conditions.  The  first,  a  Javalin  sounding  rocket,  was  launched 
by  NASA  from  Natal,  Brazil,  on  18  November  1973  at  2122  UT  [Kelley  et  al . . 
1976,  1979;  Costa  and  Kelley.  1978].  The  second,  a  Terrier-Malemute,  was 
launched  by  the  Defense  Nuclear  Agency  from  Kwajalein  Atoll  on  17  July 
1979  at  1233  UT  [Szuszczewicz  et  al..  1980;  Rino  et  al . .  1981;  Kelley  et 
al . ,  1982] .  We  will  show  that  the  spectrum  of  density  structure  observed 
in  the  F^-layer  valley  (160  to  200  km  altitude)  off  the  magnetic  equator 
is  consistent  with  that  expected  from  the  image  formation  process  (which 
is  scale  length  and  density  dependent)  driven  by  the  gravitational 
Rayleigh-Taylor  instability  operating  on  the  bottomside  F  layer  at  the 
equator  itself. 

1 .  Data  Presentation 

A  plasma  density  profile  from  the  downleg  of  the  Natal  rocket 
flight  is  reproduced  in  Figure  23(a)  with  an  expanded  plot  of  the  F- 
layer  valley  and  E-region  density  profiles  in  Figure  23(b).  The  plasma 
density  profile  was  extremely  structured  throughout  the  flight.  Near 
and  just  below  the  F-region  peak,  the  turbulence  has  been  interpreted  in 
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ALTITUDE  —  km 

200  180 


FIGURE  23  (Concluded) 

terms  of  the  non-linear  gravitational  Rayleigh-Taylor  instability  (see 
earlier  references).  The  low-altitude  (90  to  130  km)  structure  is  the 
result  of  classic  sporadic  E  layers  and  is  not  produced  by  the  "image" 
process.  Of  interest  to  the  present  study  are  the  quasi-sinusoidal 
fluctuations  in  the  F^-layer  valley  (170  to  200  km),  where  the  zero-order 
vertical  plasma  density  gradient  nearly  vanishes.  Because  of  the  finite 
dip  angle  at  Natal,  these  fluctuations  project  along  magnetic  field  lines 
to  the  magnetic  equator  at  heights  ranging  from  280  to  315  km.  Note  that 
this  altitude  range  corresponds  to  heights  where  the  local  equatorial  F 
region  was  observed  to  be  unstable  to  the  gravitational  Rayleigh-Taylor 
process.  The  geometry  is  illustrated  schematically  in  Figure  24. 


Another  example  from  the  Plumex  rocket  downleg  is  presented  in 
Figure  25.  Again,  quasi-sinusoidal  oscillations  were  observed  in  the 
F^-layer  valley  region  at  heights  that  map  to  the  magnetic  equator  in  the 
region  where  local  equatorial  spread-F  was  detected.  Conversely,  the 
electron  density  profiles  reported  by  Morse  et  al . .  [1977]  and  Mr.  S. 
Prakash  [private  communication,  1982],  which  were  obtained  very  near  the 
dip  equator  in  Peru  and  India,  respectively,  showed  no  such  structure  in 
the  F^-layer  valley.  Thus,  in  all,  five  equatorial  rocket  flights  have 
shown  E-region  structu'e  or  the  lack  thereof  in  a  manner  consistent  with 
low-altitude  structure  iieing  formed  by  the  image  process. 
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FIGURE  24  SCHEMATIC  DIAGRAM  OF  THE  MAGNETIC  FIELD  GEOMETRY 


FIGURE  25  PLUMEX  ROCKET  DOWNLEG 


A  series  of  spectra  corresponding  to  different  altitude  ranges 

from  Figure  23  are  presented  in  Figures  26  and  27.  The  Figure  26  spectra 

all  correspond  to  bottomside  equatorial  spread-F.  The  spectrum  from  the 

-2 

lowest  altitude  displays  a  form  that  varies  at  k  over  the  regime  from 
\  =  5000  m  to  X  =  50  m,  where  the  instrument  noise  level  was  reached. 

The  highest  altitude  spectrum  seems  to  have  a  somewhat  shallower  spectral 
form  at  the  long  wavelengths  (500  to  5000  m).  These  results  are  in  good 
agreement  with  the  Plumex  I  results  reported  by  Rino  et  al. .  [1981],  and 
Kellev  et  al.  [1982]. 

The  F^-layer  valley  spectra  in  Figure  27  are  quite  different 
from  those  in  Figure  26,  as  are  the  waveforms.  The  spectra  are  peaked 
between  2  to  3  Hz,  which  corresponds  to  wavelengths  of  700  to  1400  m. 

The  peaked  spectra  are  consistent  with  the  quasi-sinusoidal  nature  of  the 
waveforms.  We  show  in  the  next  section  that  such  a  peaked  spectra  should 
be  expected  on  the  basis  of  the  image  formation  process. 
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ELECTRON  DENSITY  —  km 
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TIME  —  Mcondi  liter  21:22:00.0 
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FIGURE  26  A  SERIES  OF  SPECTRA  CORRESPONDING  TO  DIFFERENT  ALTITUDE  RANGES 
FROM  FIGURE  25 
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FIGURE  26  (Continued! 
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FIGURE  26  (Concluded) 
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electron  density 


FIGURE  27  (Concluded) 


2.  Analysis 

The  mechanisms  whereby  plasma  irregularities  are  formed  in  the 
F  layer  ionosphere  can  be  studied  by  considering  the  continuity  equation: 


a,  F 

^-  +  7 

5n  ~ 


(n  v)  =  P  -  L 


(9) 


where  n  is  the  electron  density,  v  is  the  flow  velocity,  P  the  production 
rate,  and  L  the  loss  rate.  In  the  nighttime  F  layer  the  latter  two  terms 
can  be  ignored,  yielding 


.  F 

—  +  V 
5 1  ~ 


7nF  +  nF(v 


(10) 


Ignoring  variations  along  the  direction  of  the  magnetic  field,  the  primary 

2 

velocity  is  given  by  E  xB/B  ,  which  is  divergence  free,  and  thus 


V  •  7n 


(11) 


In  effect,  the  irregularities  are  formed  by  incompressibly 
mixing  an  existing  background  density  gradient  via  the  perturbation  elec¬ 
tric  fields  generated  in  the  instability  process.  These  processes  are 
reasonably  well  understood  and  well  documented  (see  recent  reviews  by 
Feler  and  Kelley  [1980]  and  Ossakow  [1981]).  For  our  purposes  here,  we 
need  only  describe  the  perturbation  electric  field  (or  equivalent  velocity) 
structure  in  the  unstable  F  layer  because  it  is  this  structure  that  maps 
to  lower  altitudes  and  drives  the  image  process. 

It  can  be  shown  from  nonlinear  Rayleigh-Taylor  instability 
theory  with  no  zero-order  zonal  electric  field  [Keskinen  and  sakow. 

1981),  that  the  electric  field  and  density  fluctuations  are  r  ated  by: 
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(12) 


where  g  is  gravitational  acceleration,  B  is  the  magnetic  field  strength, 

F  f6n"lF 

is  the  ion-neutral  collision  frequency,  and  |~ I  are  density  fluc¬ 
tuations  like  those  whose  spectra  are  shown  in  Figure  26.  Thus,  the 
velocity  turbulence  driving  image  formation  can  be  easily  specified  in 
terms  of  our  measurements  of  the  density  fluctuations.  Our  goal  then  is 
to  show  that  the  spectrum  of  image  irregularities  observed  is  the  result 
of  a  balance  between  this  scale-size-dependent  driver  and  the  image  for¬ 
mation  and  diffusive  damping  processes  that  are  also  scale-size  dependent. 
Equation  (12)  also  applies  to  midlatitude  barium  clouds  if  |g/v^j  Is 
replaced  by  (E//B),  where  E;  is  the  electric  field  measured  in  the 
neutral  frame  of  reference. 

Note  that  the  Rayleigh-Taylor  process  will  not  create  structure 
in  the  F-layer  valley  if  the  zero-order  density  gradient  vanishes.  The 
electric  field  would  still  map  to  the  region  where  density  structure  is 
observed,  but  with  no  density  gradient  the  stirring  will  not  produce 
structured  density. 

To  understand  the  formation  of  images,  we  must  examine  a  more 
complicated  version  of  the  continuity  equation.  As  will  be  discussed  in 
more  detail  below,  images  form  at  altitudes  where  the  incompressibility 
assumption  used  in  the  previous  discussion  breaks  down.  Compressibility 
arises  from  nonnegligible  Pedersen  drifts  |vp  *  —  e|  that  must  be  added 
to  the  E  x  B  flow  considered  above.  Here  ?p  is  the  Pedersen  conductivity, 
which  for  E-  and  F^-layer  heights  is  given  by 
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where  q  is  the  electron  charge,  m.  and  Q.  are  the  ion  mass  and  gyrofre- 

E  1  1 

quency,  respectively,  «»ud  is  the  ion  neutral  collision  frequency  at 
the  image  altitude.  Throughout  what  follows,  the  superscript  E  will 
denote  quantities  at  the  image  altitude  even  though  we  are  concentrating 
principally  on  the  F^-layer  valley  region.  Additional  complications  to 
the  continuity  equation  include  the  effects  of  relatively  rapid  cross¬ 
field  plasma  diffusion  and  chemical  recombination.  If  we  include  these 
effects,  Eq.  (10)  becomes 
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where  we  have  ignored  any  horizontal  variation  in  electron  density  but 
have,  for  completeness,  included  a  vertical  gradient  term  that  is  weak  in 

the  valley  region.  D  is  the  cross-field  plasma  diffusion  rate,  and  a  is 

x  E 

the  effective  recombination  coefficient.  Dividing  through  by  nQ,  we 

obta;n  an  expression  for  the  linear  growth  rate  of  the  image  striation, 
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[15(a)] 


Thus,  growth  of  images  will  occur  if 


V  »  V  4*  V  ■  V  •  Y  0 

image  stir  Tconv  diff  chem 


[15(b)] 


Before  discussing  the  physical  significance  of  each  of  the  terms, 
we  should  recall  that  the  electric  field,  E,  is  the  perturbation  electric 
field  associated  with  the  driving  instability  operating  on  the  bottomside 
F  layer  and  mapped  to  the  lower  altitude  image  region.  For  the  electric 
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field  fluctuations  to  map,  they  must  be  of  scale  sizes  larger  than 
.several  hundred  meters  [Farlev.  1959,  I960].  Also,  for  the  steady  state 
conditions  considered  here,  images  must  form  at  the  same  scale  size  as 
the  driving  F-layer  instability  striations  [see  Vickrey  and  Kelley.  1982]. 
An  additional  subtlety  of  Eqs.  (13)  and  [15(a)]  is  that  we  have  assumed 
in  the  diffusion  term,  to  be  independent  of  density  so  that  we 

could  bring  it  outside  the  divergence  operator  In  reality,  D  depends 
on  electron  density  through  electron-ion  collisions,  but  at  the  altitudes 
of  present  interest  this  dependence  is  negligible  compared  with  the 
dependence  on  ion-neutral  collisions. 

The  physical  processes  that  compete  to  determine  whether  an 
image  forms  at  a  given  scale  size  and  background  density  may  be  described 
as  follows.  The  stirring  term  yst£r  i-n  Eq.  [15(b)]  can  produce  an  image 
if  plasma  of  different  densities  (from  different  heights  in  the  equator¬ 
ial  situation)  is  mixed  on  a  time  scale  that  is  short  compared  with  the 
time  it  takes  for  the  plasma  to  recombine  or  diffuse  away.  This  time 
scale  is  determined  by  the  magnitude  of  the  mixing  velocity  and  the 
steepness  of  the  density  gradient.  In  the  present  situation,  the  density 
gradient  of  interest  is  the  vertical  scale  height  at  the  image  altitude. 
However,  because  this  scale  height  Ly  was  so  large  in  the  present  case 
[L'v  ~  60  km  in  Figure  22(a)],  this  term  is  important  only  at  very  large 
scales.  It  may  be  noteworthy  for  other  situations,  however,  that  the 
scale-size  dependence  of  Ygt£r  Is  the  same  as  that  of  the  driving  velo¬ 
city  fluctuations  in  the  unstable  F  layer. 

The  second  term  in  Eq.  [15(b)],  Yconv»  represents  the  formation 
of  an  image  caused  by  a  converging  velocity  field  operating  on  the  com¬ 
pressible  ion  gas.  The  electron  gas  is  still  incompressible  at  image 
altitudes,  and  thus  electrons  must  flow  up  and  down  along  magnetic  field 
lines  to  maintain  charge  balance  [see  discussion  by  Vickrev  and  Kelley. 
1982a] .  Note  that  the  divergence  operator  introduces  an  additional 
scale-size  dependence  to  the  v  term.  Small-scale  fluctuations  in  the 
velocity  field  form  images  more  readily  than  do  large-scale  fluctuations 
of  the  same  magnitude  [Vickrey  and  Kelley.  1982a] .  As  was  the  case  with 
the  stirring  term,  however,  the  convergence  of  plasma  must  proceed  at  a 
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rate  faster  than  its  recombination  rate  in  order  for  an  image  to  form. 
Thus,  as  the  background  density  is  increased,  smaller  scale  image  growth 
will  be  favored. 


The  term  vcjiem  *-n  Eq.  [15(b)]  represents  chemical  recombination 
of  the  image  region  plasma.  At  altitudes  where  the  ion  gas  is  compress¬ 
ible,  the  ionospheric  ion  composition  tends  to  be  dominated  by  molecular 

ions  (i.e.,  N0+  and  ot)..  The  recombination  rate  for  these  species  varies 
2  ^ 

as  cmE  ,  which  produces  a  strong  density  dependence  in  contrast  to  the  F- 

layer  plasma  where  the  driving  instability  operates.  At  F-layer  heights 

+  p  2 

the  recombination  rate  of  0  is  proportional  to  n  instead  of  nF  . 

1  ^nF 

Hence,  the  time  constant  for  recombination,  "TT",  is  independent  of 

n0  ot 

electron  density. 


Now  let  us  examine  the  scale-length  dependence  of  the  image 
growth  rate  expression.  Making  use  of  Eq.  (12)  we  can  Fourier  Transform 
Equation  [15(a)]  to  obtain 
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where  P(k)  =*  — —  (k).  The  spectra  in  Figure  25  are  well  represented  by 
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corresponding  to  k£ 


6.2  x  10“3m_1. 
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Figure  28  shows  plots  of  V ima^(k)  versus  k  for  typical  F-layer 

valley  conditions  corresponding  to  the  observations  discussed  above,  i.e., 

vE  /n.  =  0.1,  g/vF  =  100  m  s'1,  D  =  10  m2  s'1,  and  a  =  5  x  10'8cm3S-1 
in  i  in  i 

[Vickrey  et  al . .  1982].  It  is  important  to  point  out  that  the  detailed 
shape  of  the  curve  and  the  magnitude  of  y(k)  depend  on  the  choice  of  the 
integration  interval  Ak  in  Eq.  (16).  To  produce  the  curves  we  have 
chosen  values  of  Ak  that  are  roughly  equal  to  the  bandwidth  of  the  observed 
irregularity  spectrum.  This  interval  represents  a  "bandpass"  of  driving 
fluctuation  wavenumbers  about  the  image  scale  size  of  interest.  Because 
the  choice  of  integration  interval  is  arbitrary,  one  must  not  attach  too 
much  significance  to  the  magnitude  of  the  growth  rates  in  Figure  28. 
Nonetheless,  the  theoretical  prediction  of  a  narrow  band  of  preferred 
image  scale  size  is  very  significant.  A  comparison  of  the  curves  labeled 
(a)  and  (b)  in  Figure  28  illustrates  the  effect  of  doubling  the  integra¬ 
tion  interval  Ak.  As  one  would  expect,  the  growth  rate  increases  for  all 
scales  with  the  integration  interval  Also,  the  most  preferred  value  of 
k  increases,  but  only  by  about  20  percent.  Thus,  the  prediction  of  a 
preferred  scale  is  not  highly  sensitive  to  the  integration  interval  chosen. 

The  scale- size  dependence  of  the  image  growth  rate  also  depends, 

2 

of  course,  on  the  irregularity  spectrum  of  bottomside  spread-F,  P  (k), 
which  is  assumed  to  drive  the  image  process.  To  examine  the  sensitivity 
of  Yj^aggCk)  to  the  assumed  spectral  shape,  we  have  halved  and  doubled 
the  value  of  k£  in  Eq.  (17)  to  construct  curves  (c)  and  (d)  in  Figure  28. 
The  integration  interval  was  the  same  as  that  for  curve  (a).  As  can  be 
seen  in  the  figure,  even  these  rather  dramatic  changes  in  the  driving 
spectral  characteristics  do  not  drastically  change  the  preferred  scale 
size  produced.  Our  best  estimate  of  the  conditions  appropriate  to  the 
observations  presented  above  corresponds  to  curve  (a),  which  predicts  a 
preferred  image  scale  of  ~  750  m.  This  prediction  is  in  good  agreement 
with  the  preferred  irregularity  scale  sizes  shown  in  Figure  27.  One 
must  keep  in  mind,  however,  that  although  a  peaked  growth  rate  is  concep¬ 
tually  pleasing  in  a  process  that  ultimately  results  in  a  narrow  spectrum 
of  density  irregularities,  there  is  no  fundamental  reason  why  the  fully 
developed  (nonlinear)  density  spectrum  should  necessarily  match  the 
linear  growth  rate  exactly. 
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WAVENUMBER  —  m'1 

FIGURE  28  WAVELENGTH  DEPENDENCE  OF  THE  IMAGE  GROWTH  RATE  FOR  FrREGlON 
VALLEY  ALTITUDES.  The  various  curves  illustrate  the  sensitivity  of  the  results  to 
changes  in  the  integration  interval  A/k  and  driving  spectrum  P(k). 


To  determine  whether  the  magnitude  of  the  predicted  image 
growth  rates  is  significant,  we  can  compare  it  with  the  linear  growth 
rate  of  the  driving  process,  YR-p- 


v 

conv 

yRT 


(18) 


where  Lp  is  the  zero-order  vertical  gradient  on  the  bottomside  of  the 
F  layer,  and 
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rk  +Ak/2 

f°  2 
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(19) 


Since  (— )  Ak  is  about  5  percent  in  the  bandpass  interval  of 
£ 

the  observed  images,  (v.  /D.)  ~  0.1,  L„  ~  20  km,  and  the  value  of  k  at 

ini  F_2-l  0 

the  peak  of  the  image  growth  rate  is  ~  10  m  ,  then  R  »  1.  Thus,  the 

growth  of  images  should  occur  on  a  similar  time  scale  to  that  of  the 

"parent"  striations.  We  show  below  that  v.  is  even  faster  at  lower 

r  image 

altitudes  for  comparable  electron  densities.  Thus,  since  the  occurrence 
of  bottomside  spread-F  is  well  known  and  documented,  the  formation  of 
image  striations  should  not  be  too  surprising. 


Although  we  have  no  direct  observations  of  image  formation  in 
the  E  region,  the  process  should  proceed  more  rapidly  there  (provided 
the  density  is  low  enough,  e.g.,  Vickrey  and  Kelley.  1982a)  because  the 
ion  gas  is  more  compressible  at  E-region  heights  than  at  F-region  heights. 
In  Figure  29  we  have  calculated  the  image  growth  rate  V.  (k)  for  E- 

2  -1  _7  .3 

region  conditions;  i.e.,  v^/Q  *  1,  D  ■  100  m  s  ,  and  a  *>  2  x  10  cm  5 
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FIGURE  29  WAVELENGTH  DEPENDENCE  OF  THE  IMAGE  GROWTH  RATE  FOR  E-REGION 
ALTITUDES 


0 


8 


[Vickrey  et  al. ,  1982] .  The  integration  interval  and  kc  are  the  same  as 
Figure  28  curve  (a).  We  find,  as  at  higher  altitudes,  a  preferred  image 
scale  size  of  ~  700  m  is  produced.  In  this  case,  however,  the  image 
growth  rate  is  approximately  two  orders  of  magnitude  faster  than  at  F  - 
layer  altitudes  provided  the  electron  density  is  comparable. 

We  have  presented  here  what  we  believe  to  be  the  first  direct 
evidence  for  the  formatidri  of  image  irregularities  in  the  ionospheric 
medium.  We  have  also  computed  a  linear  growth  rate  for  the  structures 
and  find  that  the  k  dependence  matches  the  observed  spectrum  quite  well. 

The  case  we  have  studied  involves  images  of  F-region  irregular¬ 
ities  created  by  the  Rayleigh-Taylor  process  in  equatorial  spread-F.  A 
similar  analysis  should  apply,  however,  to  midlatitude  barium  cloud 
striations  that  should  also  create  E-region  images  in  the  same  scale-size 
regime  as  presented  here.  For  example,  Kelley  et  al.  [1979]  showed  that 
power  spectral  measurements  in  a  striated  barium  cloud  matched  quite 
closely  the  bottomside  spread-F  spectrum  measured  in  the  Natal  flight 
discussed  here  in  both  absolute  magnitude  and  scale  size.  Thus,  E-region 
structures  should  have  formed  with  comparable  scales  to  those  reported 
here  and  with  growth  rates  on  the  order  of  seconds. 

At  auroral  latitudes  the  application  is  not  so  clear  because 
the  chemical  damping  term,  crao,  is  larger  than  at  midlatitudes  owing  to 
particle  precipitation.  Exceptions  might  be  in  the  winter  polar  cap  or 
in  the  region  just  equatorward  of  the  diffuse  auroral  zone  where  the  E- 
region  density  should  be  comparable  to  midlatitude  values. 


IV  SUGGESTIONS  FOR  FUTURE  WORK 


The  results  we  have  outlined  from  the  past  year's  work  have  spanned 
a  large  area  of  theoretical  and  experimental  physics.  It  is  clear  that 
the  area  least  understood  in  both  the  natural  and  the  nuclear  environ¬ 
ments  is  high-latitude  plasma  structure  phenomenology.  During  this  year 
we  intend  to  concentrate  principally  on  this  area,  making  heavy  use  of 
data  from  the  HILAT  satellite  in  conjunction  with  incoherent-scatter 
radar  and  spaced-receiver  measurements.  We  will  develop  theoretical 
models  to  guide  the  interpretation  of  the  experimental  data;  where  poss¬ 
ible,  we  will  use  these  models  to  simulate  nuclear  effects  at  high  lati¬ 
tudes.  One  simple  but  nevertheless  fundamental  question  that  we  will 
address  is  the  role  of  the  magnetosphere  in  prolonging  or  enhancing  the 
development  of  structure  in  the  nuclear-disturbed  environment.  Such 
prolongation  might  be  expected  because  destabilizing  magnetospheric  field- 
aligned  currents  will  be  intensified  at  the  edges  of  the  strong  conductivity 
enhancement  produced  by  the  detonation.  Another  direct  application  of  our 
magnetospheric  convection  model  will  be  to  predict  where  nuclear-produced 
debris  ionization  will  be  distributed  at  very  late  times.  The  redistribu¬ 
tion  will,  of  course,  be  a  strong  function  of  the  latitude  and  local  time 
of  the  detonation. 

Further  refinement  and  verification  of  our  developing  models,  as 
well  as  those  of  others,  will  be  carried  out  through  extensive  analysis 
of  data  on  naturally  occurring  F-layer  structure  from  HILAT.  We  have 
previously  presented  a  framework  for  understanding  the  morphology  of 
naturally  occurring  high- latitude  irregularities  [Vickrey.  1982;  DNA  La 
Jolla  Conference,  1982;  Vickrev  and  Kellev.  1982;  Kellev  et  al. .  1982]. 

The  framework  consists  of  accurately  specifying:  the  latitude-local¬ 
time  "source  function"  of  irregularity  production,  the  (scale-size  depen¬ 
dent)  lifetime  of  irregularities  once  produced,  and  the  redistribution  of 
structured  plasma  during  its  lifetime  by  magnetospheric  convection. 

Although  this  framework  is  straightforward,  each  individual  element  is 
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very  complex  but  can  be  addressed  using  HILAT  data.  Listed  below  are 
brief  descriptions  of  some  of  the  studies  we  intend  to  begin  in  collabora¬ 
tion  with  other  DNA  researchers  during  the  next  year.  This  list  is  grouped 
according  to  our  irregularity  morphology  framework.  Obviously,  the  list 
is  too  extensive  to  be  completely  fulfilled  within  one  year,  but  we  hope 
to  make  substantive  progress  in  most  of  these  foreseeable  problem  areas. 

A.  Irregularity  Production  Studies 

1 .  Global  Precipitation  Source  Function 

It  is  important  to  characterize  accurately  the  latitude-local- 
time  variations  of  precipitating  particle  flux  as  a  function  of  particle 
energy.  The  changes  in  precipitation  characteristics  with  magnetic 
activity,  in  the  interplanetary  magnetic  field  direction,  in  solar  wind 
speed,  and  in  the  general  state  of  the  magnetosphere  should  be  specified. 
This  study  can  be  accomplished  in  large  part  through  routine  analysis  and 
careful  sorting  of  the  particle  detector  data  from  HILAT.  However,  early 
in  the  satellite  mission  the  2-D  vacuum  ultra  violet  (VUV  )  images  should 
be  calibrated  against  the  particle  detector  and  incoherent- scatter 
measurements  so  that  they  can  be  interpreted  both  in  terms  of  the  source 
function  itself  and  perhaps  in  specifying  the  state  of  the  magnetosphere. 

2.  Spatial  Structure  of  Precipitation 

To  date,  the  spectral  characteristics  of  particle  precipitation 
have  not  been  examined  in  detail — particularly  at  the  soft  energies  res¬ 
ponsible  for  F-region  ionization.  It  is  not  known,  therefore,  if  struc¬ 
tured  precipitation  can  directly  produce  structured  F-layer  ionization  at 
kilometer  scales  and  smaller,  or  if  precipitation  is  only  responsible  for 
the  large-scale  "seed"  features  on  which  instabilities  operate.  To 
examine  this  question,  the  same  spectral  analysis  techniques  that  will  be 
applied  to  the  in-situ  density  data  (as  well  as  the  beacon,  driftmeter, 
and  magnetometer  data-- see  below)  should  be  applied  to  the  particle 
detector  data.  In  this  way  the  direct  contribution  of  the  spatial  preci¬ 
pitation  spectrum  to  the  in-situ  density  irregularity  spectrum  can  be 
immediately  assessed. 
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3. 


Global  Conductivity  Distribution 


The  particle  detector  and,  to  some  extent,  the  VUV  imager  can 
be  analyzed  (as  in  the  first  study  above)  to  give  the  latitude-local-time 
distribution  of  precipitating  particle  flux  as  a  function  of  energy. 

Using  a  model  neutral  atmosphere,  this  information  can  be  used  to  derive 
the  ionospheric  E-region  density  profile,  and,  hence,  the  conductivity,  £ . 

The  importance  of  determining  the  E-region  conductivity  is  that  it 
strongly  affects  both  the  growth  rate  of  convective  instabilities  and  the 
diffusive  decay  of  irregularities  once  produced. 

4.  Global  Characterization  of  Density  Irregularity  Spectrum 

Some  irregularity  sources  and  sinks  may  be  highly  scale-size 
dependent.  To  examine  the  interaction  of  these  mechanisms,  the  entire 
spectrum  of  irregularities  must  be  studied  as  a  whole.  Moreover,  non¬ 
linear  plasma  processes  can  couple  energy  between  different  scale-size 
regimes.  Therefore,  a  careful  spectral  analysis  of  the  in-situ  density 
data  is  essential  so  that  the  latitude-local-time  variations  in  the  irregu¬ 
larity  spectrum,  S(k),  can  be  characterized.  Evaluation  of  propagation 
effects  on  communications  and  surveillance  systems  begins  with  this 
characterization. 

5 .  Convection  Velocity  Shears 

Structured  electric  fields  can  produce  structured  plasma  density 
in  two  principal  ways.  First,  where  the  electric  field,  E,  is  large,  the 
enhanced  ion  temperature  due  to  Joule  heating  affects  the  F-region  chemistry 
in  such  a  way  as  to  produce  a  density  depletion.  At  present,  we  do  not 
know  over  what  scale-size  regime  this  process  might  operate.  Second,  a 
structured  electric  field  can  produce  structured  density  by  simply  mixing 
flux  tubes  of  different  densities.  Thus,  a  turbulent  magnetospheric 
electric  field  can  produce  much  the  same  result  as  a  local  convective 
instability.  Therefore,  it  is  important  that  the  magnitude,  frequency  of 
occurrence,  and  global  distribution  of  velocity  shears  be  characterized 
and  correlated  with  regions  of  density  structure. 
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Relationships  Between  E  and  VN 


An  examination  of  the  relationship  between  E  and  the  density 
gradient  (VN)  should  give  an  indication  of  whether  structure  at  a  given 
scale  size  is  being  driven  by  an  instability  or  is  diffusing.  In  the 
latter  case,  theoretical  predictions  of  the  relationship  between  6E  and 
6N  under  varying  E-region  conductivity  conditions  can  be  tested;  in  the 
instability  case,  the  relationship  between  these  quantities  can  help 
identify  wave  modes;  e.g.,  fluted  or  nonfluted. 

7.  Global  and  Small-Scale  Parallel  Current  Patterns 

The  field-aligned  current  magnitude  ( J u ) ,  is  the  driving  term 
in  the  growth  rate  expression  for  the  current-convective  instability. 

The  sense  of  the  large-scale  global  patterns  of  field-aligned  current 
has  been  studied  extensively  with  TRIAD  data,  but  the  magnitude  of  the 
currents  has  not  been  studied.  Moreover,  interest  has  been  focused  in 
the  past  on  the  global  morphology  of  current  systems.  Although  these 
large-scale  features  provide  a  permanent  destabilizing  influence,  small- 
scale  field-aligned  currents  may  be  much  more  intense.  It  thus  seems 
worthwhile  to  apply  the  same  spectral  analysis  technique  used  on  the 
density,  velocity,  and  precipitating  particle  data  to  the  magnetometer 
data.  In  that  way  the  magnitude  and  intensity  of  small-scale  currents 
can  be  assessed  statistically. 

8.  Global  Empirical  Modeling  of  E.  5.  J u .  and  U 

Indications  are  that  the  pertinent  instability  growth  rates  in 
the  polar  regions  are  faster  than  the  classical  diffusive  lifetime  of 
irregularities.  This  means  that  growth  of  structure  on  a  given  set  of 
flux  tubes  is  cumulative,  to  some  extent,  as  those  flux  tubes  convect 
around  the  polar  regions.  The  amount  of  structure  on  those  flux  tubes 
therefore  depends  on  their  past  history.  An  assessment  of  the  destabi¬ 
lizing  influences  suffered  during  that  time  can  be  assessed  only  by 
modeling  the  global  patterns  of  E,  S,  J u ,  and  the  neutral  wind,  U. 
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Comparison  of  TEC  and  Scintillation  Boundary 

The  extensive  Wideband  data  set  shows  that  the  equatorward  TEC 
boundary  is  consistently  associated  with  a  region  of  enhanced  scintilla¬ 
tion-producing  irregularities.  This  relationship  merits  further  investi¬ 
gation  with  the  more  extensive  diagnostics  on  board  HILAT.  It  may  be 
possible,  for  example,  to  relate  the  equatorward  TEC  boundary  directly 
to  the  equatorward  magnetospheric  convection  boundary. 

B.  Irregularity  Lifetime  Studies 

1.  Global  Conductivity  Distribution 

When  the  E-region  Pedersen  conductivity  is  high,  the  classical 
F-region  cross-field  diffusion  rate  can  be  enhanced  by  an  order  of  magni¬ 
tude.  Also,  the  pertinent  instability  growth  rates  scale  by  the  ratio 
of  F-region  to  E-region  Pedersen  conductivity.  Thus,  the  above  Global 
Conductivity  Distribution  study  is  pertinent  to  two  of  our  framework 
elements. 

2.  Global  Characterization  of  Density  Irregularity  Spectrum 

To  examine  irregularity  sources  and  sinks,  the  entire  spectrum 
of  irregularities  must  be  studied  as  a  whole,  using  data  from  both  the 
drift  meter  and  multifrequency  phase  propagation  observations.  This  may 
reveal  scale-size  dependent  drivers  and  damping  mechanisms.  There  are 
some  indications  that  anomalous  diffusion  operates  at  high  latitudes,  but 
this  assumption  is  based  on  examining  irregularities  at  only  one  scale 
size.  If  the  diffusion  is  due  to,  say,  drift  waves,  then  it  may  be  highly 
scale-size  dependent. 

3.  Relationships  Between  E  and  VN 

As  in  the  source  studies  just  mentioned,  if  anomalous  diffusion 
is  caused  by  waves  in  the  plasma,  then  the  relationship  between  6E  and  6N 
will  help  to  identify  those  waves.  If  anomalous  diffusion  is  driven  by 
microinstabilities  (which  produce  waves  with  perturbed  E  on  the  scale 
size  of  the  ion  gyroradius),  then  perhaps  rocket  measurements,  in  addition 

to  HILAT  measurements,  would  be  required  to  resolve  the  waves. 
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Global  Empirical  Modeling  of  E.  £.  Jn .  and  U 

As  discussed  in  A. 7  above,  structured  growth  is  cumulative.  By 
|  the  same  token,  if  flux  tubes  convect  into  regions  where  the  configuration 

I  of  E,  U,  and  VN  is  stabilizing,  then  structure  should  be  damped. 

C.  Convection-Driven  Redistribution  of  Structured  Plasma 

[ 

The  high-latitude  convection  pattern  is  known  to  depend  strongly  on 
the  orientation  of  the  interplanetary  magnetic  field,  solar  wind  speed, 
and  so  on.  This  behavior  needs  further  characterization  to  ultimately 
achieve  predictive  capability.  Not  only  does  convection  redistribute 
structured  plasma  (and  perhaps  play  a  role  in  determining  irregularity 
geometry),  but  slippage  between  the  convective  flow  and  the  neutral  wind 
can  produce  the  structure.  > 

An  additional  and  needed  system  specification,  as  well  as  a  helpful 
clue  for  understanding  nonlinear  structure  evolution,  is  the  geometry  of 
irregularities.  This  geometry  is  undoubtedly  scale-size  dependent  and 
cannot  be  determined  from  the  satellite  in-situ  data  alone.  Its  deter¬ 
mination  will  require  different  measurements  for  different  scale  sizes; 
e.g.,  incoherent  scatter  measurements  for  \  >  10  km,  and  spaced-receiver 
measurements  for  \  <  1  km. 

Spaced-receiver  measurements  from  HILAT  will  extend  the  coverage  of 
kilometer-scale  irregularity  anisotropy  measurements  to  other  local  times 
that  Wideband  did  not  sample.  Extension  of  the  latitudinal  coverage  over 
that  of  Wideband  is  also  possible.  The  comparison  of  data  taken  from 
HILAT  with  data  currently  being  analyzed  from  quasi-stationary  satellites 
will  also  help  elucidate  the  differences  between  temporal  and  spatial  de- 
correlation  of  scintillation-producing  structure. 

Recent  measurements  (as  well  as  simple  theory)  [Vickrey,  1982]  suggest 
that  irregularity  anisotropy  may  be  strongly  influenced  by  the  magneto- 
spheric  convection  pattern.  Coordinated  measurements  taken  by  spaced- 
receivers  at  Sondrestrom,  the  ion  drift  meter  aboard  HILAT,  and  incoherent- 
scatter  radar  can  examine  this  relationship  and  its  scale-size  dependence. 
We  also  intend  to  conduct  experiments  in  conjunction  with  tb®  EISCAT 
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incoherent-scatter  facility.  Because  of  its  unique  capabilities  (in 
particular,  its  tristatic  configuration.  Barker  coding,  lower  frequency 
and  increased  sensitivity,  compared  with  those  of  the  Chatanika  facility), 

E IS CAT  may  be  able  to  resolve  structure  in  the  scale- size  regime  not 
covered  by  the  Sondrestrom  radar/ spaced- receiver  measurements;  i.e.,  1  to 
10  km.  The  tristatic  configuration  may  also  allow  very  direct  comparison 
of  irregularity  anisotropy  and  convection  streamlines.  Moreover,  because 
the  density  configuration  and  velocity  field  can  be  measured  simultaneously 
with  the  same  resolution,  one  can  immediately  assess  the  growth  rate  of 
structure  for  comparison  to  theory. 
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Ooc  Control 

DEPARTMENT  OF  THE  NAVY 

Nival  Ocean  System  Center 
ATTN:  Code  532 

ATTN:  Code  5323,  J.  Ferguson 

ATTN:  Code  5322,  N.  Paulson 

Naval  Research  Lab 

ATTN:  Code  4720,  0.  Davis 

ATTN:  Code  7500.  8.  Maid 

ATTN:  Code  4780,  S.  Ossakow 

ATTN:  Code  4700 

ATTN:  Code  4780 

ATTN:  Code  6700 

ATTN:  Code  4187 

ATTN:  Code  7950,  J.  Goodman 

Office  of  Naval  Research 

ATTN:  Code  414,  G.  Joiner 
ATTN:  Code  412,  W.  Condell 

DEPARTMENT  OF  THE  AIR  FORCE 

Air  Force  Weapons  Lab 
ATTN:  SUL 

ATTN:  NTYC 

ATTN:  NTN 


DEPARTMENT  OF  THE  AIR  FORCE  (Continued) 

Air  Force  Geophysics  Lab 

ATTN:  OPR,  H.  Gardiner 
ATTN:  OPR-1 
ATTN:  LKB,  K.  Champion 
ATTN:  CA,  A  Stafr 
ATTN:  PHY,  J.  Buchau 
ATTN:  R.  Babcock 
ATTN:  R.  O'Neil 

Air  Force  Tech  Applications  Ctr 
ATTN:  TN 

Air  Force  Wright  Aeronautical  Lab 
ATTN:  W.  Hunt 
ATTN:  A.  Johnson 

Rom  Air  Oevelopaient  Center 
ATTN:  OCS,  V.  Coyne 
ATTN:  TSLD 

Rose  Air  Development  Center 
ATTN:  EEP,  J.  Rasmussen 

OTHER  GOVERMtNT  AGENCIES 

Department  of  Commerce 
ATTN:  R.  Grubb 

Institute  for  Telecoami  Sciences 
ATTN:  L.  Berry 
ATTN:  A.  Jean 
ATTN:  W.  Utlaut 

DEPARTMENT  OF  ENER6Y  CONTRACTORS 

Los  Alamos  National  Lab 

ATTN:  MS  670.  J.  Hopkins 
ATTN:  P.  Keaton 
ATTN:  MS  664,  J.  Zlnn 
ATTN:  T.  Kunkle,  ESS-5 
ATTN:  R.  Jeffries 
ATTN:  D.  Simons 
ATTN:  J.  Wolcott 

Sandla  National  Labs 

ATTN:  D.  Thombrough 
ATTN:  Tech  Lib  3141 
ATTN:  D.  Dahlgren 
ATTN:  Space  Project  Dlv 
ATTN:  Org  4231,  T.  Wright 
ATTN:  Org  1250,  W.  Broun 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corp 

ATTN:  R.  Slaughter 
ATTN:  J.  Straus 
ATTN:  V.  Josephson 
ATTN:  T.  Salmi 
ATTN:  D.  Olsen 
ATTN:  K.  Cho 
ATTN:  I.  Garfunkel 
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DEPARTMENT  Of  OEFENSE  CONTRACTORS  (Continued) 

SOM  Corp 

ATTN:  L.  Jacobs 
ATTN:  T.  Neighbors 

Berkeley  Research  Associates,  Inc 
ATTN:  C.  Prattle 
ATTN:  J.  Workman 
ATTN:  S.  Brecht 

Charles  Stark  Draper  Lab,  Inc 
ATTN:  D.  Cox 
ATTN:  A.  Tetewskl 
ATTN:  J.  Gilmore 

Cornell  University 

ATTN:  D.  Farley  Jr 
ATTN:  M.  Kelly 

EOS  Technologies,  Inc 
ATTN:  B.  Gabbard 

General  Research  Corp 
ATTN:  8.  Bennett 

HSS,  Inc 


ATTN: 

D.  Hansen 

Institute  for  Defense  Analyses 

ATTN: 

J.  Aeln 

ATTN: 

E.  Bauer 

ATTN: 

H.  Wolfhard 

ATTN: 

H.  Gates 

JAYCOR 

ATTN: 

J.  Sperling 

Kanan  Tempo 

ATTN: 

0AS1AC 

ATTN: 

B.  Gembltl 

ATTN: 

J.  Devore 

ATTN: 

K.  Schwartz 

ATTN: 

W.  McNamara 

ATTN: 

J.  Devore 

Kaamn  Tempo 

ATTN: 

OASIAC 

M.I.T.  Lincoln  Lab 

ATTN: 

0.  Towle 

Mission  Research  Corp 

ATTN: 

R.  Hendrick 

ATTN: 

F.  Fajen 

ATTN: 

R.  Kllb 

ATTN: 

R.  Bogusch 

ATTN: 

C.  Lauer 

ATTN: 

S.  Gutsche 

ATTN: 

D.  Knapp 

ATTN: 

R.  Blgonl 

ATTN: 

Tech  Library 

ATTN: 

G.  NcCartor 

ATTN: 

F.  Gutgllano 

DEPARTMENT  Of  DEFENSE  CONTRACTORS  (Continued) 

Pacific-Sierra  Research  Corp 
ATTN:  E.  Field  Jr 
ATTN:  H.  Brode,  Chairman  SAGE 
ATTN:  F.  Thomas 

Photometries,  Inc 

ATTN:  I.  Kofsky 

Physical  Dynamics,  Inc 
ATTN:  E.  Fremouw 

Physical  Research,  Inc 
ATTN:  R.  Dellberls 

RAD  Associates 

ATTN:  R.  Letevler 
ATTN:  U.  Karzas 
ATTN:  W.  Wright 
ATTN:  R.  Turco 
ATTN:  M.  Gantsweg 
ATTN:  H.  Ory 
ATTN:  C.  Grelflnger 
ATTN:  F.  Gilmore 

RAD  Associates 


ATTN:  B. 

Yoon 

SRI  International 

ATTN:  V. 

Gonzales 

ATTN:  U. 

Chesnut 

ATTN:  D. 

McDaniels 

ATTN:  M. 

Baron 

ATTN:  R. 

Leadabrand 

ATTN:  G. 

Price 

ATTN:  0. 

Nell son 

ATTN:  A. 

Burns 

ATTN:  W. 

Jaye 

ATTN:  J. 

Petrlckes 

ATTN:  G. 

Smith 

4 

cy 

ATTN:  R. 

Tsunoda 

4 

cy 

ATTN:  J. 

Vickrey 

4 

cy 

ATTN:  R. 

Livingston 

4 

cy 

ATTN:  C. 

Rlno 

4 

cy 

ATTN:  R. 

Robinson 

4 

cy 

ATTN:  M. 

Cousins 

4 

cy 

ATTN:  N. 

Walker 

Stewart  Radiance  Lab 
ATTN:  J.  Ulwlch 


Technology  International  Corp 
ATTN:  W.  Boqulst 

Toyon  Research  Corporation 
ATTN:  J.  Isle 
ATTN:  J.  Garbarlno 

Vlsldyne,  Inc 

ATTN:  0.  Shepard 
ATTN:  W.  Reldy 
ATTN:  C.  Humphrey 
ATTN:  J.  Carpenter 
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